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We generalize to the case of spinning black holes a recently introduced "eflective one-body" 
approach to the general relativistic dynamics of binary systems. We show how to approximately 
map the conservative part of the third post-Newtonian (3 PN) dynamics of two spinning black 
holes of masses mi, m2 and spins S\, S2 onto the dynamics of a non-spinning particle of mass 
= mi m2/(mi -|- 7712) in a certain effective metric g'^{x^; M, v, a) which can be viewed either as 
a spin-deformation (with deformation parameter a = Ses/M) of the recently constructed 3 PN 
effective metric g'^^{x^; M, v), or as a j^-deformation (with comparable-mass deformation parameter 
V = mim2/(mi + 7712)^) of a Kerr metric of mass M = mi -I- m2 and (effective) spin Sc« = 
(1 -I- 3m2/(4mi)) Si + (1 + 3mi/(4m2)) 82- The combination of the effective one-body approach, 
and of a Pade definition of the crucial effective radial functions, is shown to define a dynamics with 
much improved post-Newtonian convergence properties, even for black hole separations of the order 
of 6 GM/c? . The complete (conservative) phase-space evolution equations of binary spinning black 
hole systems are written down and their exact and approximate first integrals are discussed. This 
leads to the approximate existence of a two-parameter family of "spherical orbits" (with constant 
radius), and of a corresponding one-parameter family of "last stable spherical orbits" (LSSO). These 
orbits are of special interest for forthcoming LIGO/VIRGO/GEO gravitational wave observations. 
The binding energy and total angular momentum of LSSO's are studied in some detail. It is argued 
that for most (but not all) of the parameter space of two spinning holes the approximate (leading- 
order) effective one-body approach introduced here gives a reliable analytical tool for describing 
the dynamics of the last orbits before coalescence. This tool predicts, in a quantitative way, how 
certain spin orientations increase the binding energy of the LSSO. This leads to a detection bias, in 
LIGO/VIRGO/GEO observations, favouring spinning black hole systems, and makes it urgent to 
complete the conservative effective one-body dynamics given here by adding (resummed) radiation 
reaction effects, and by constructing gravitational waveform templates that include spin effects. 
Finally, our approach predicts that the spin of the final hole formed by the coalescence of two 
arbitrarily spinning holes never approaches extremality. 



I. INTRODUCTION 



The most promising candidate sources for the LIGO/VIRGO/GEO/. . . network of ground based gravitational wave 
(GW) interferometric detectors are coalescing binary systems made of massive (stellar) black holes [0-^- Signal to 
noise ratio (SNR) estimates suggest that the first detections will concern black hole binaries of total mass > 25Mq. 
Modelling the GW signal emitted by such systems poses a difficult theoretical problem because the observationally 
most "useful" part of the gravitational waveform is emitted in the last ~ 5 orbits of the inspiral, and during the 
"plunge" taking place after crossing the last stable circular orbit. The transition between (adiabatic) inspiral and 
plunge takes place in a regime where the two bodies are moving at relativistic speeds {v/c ^ 1/V6 ~ 0.4) and where 
their gravitational interaction becomes (nearly by definition) highly non-linear {GM/c^r ^ 1/6). 

Several authors (notably |^,D) have taken the view that the modelling of this crucial transition between inspiral 
and plunge is (in the general case of cornparable-mass systems) beyond the reach of analytical tools and can only 
be tackled by (possibly special-purpose |^) numerical simulations. By contrast, other authors |]7|-p|,p|, p^|jTl| ] have 
introduced new "resummation methods" to improve the analytical description of the last few GW cycles near this 
transition and have argued that these resummed analytical results gave a reliable description of the gravitational 
physics near the transition. The purpose of the present paper is to further the latter approach by generalizing the 
(resummed) "effective one-body" (EOB) methods introduced in pj-p^ to the case of binary systems of spinning black 
holes. Before doing this, we wish to clarify what is the rationale for arguing that the "resummed" analytical approach 
can describe the last stages of inspiral and the transition between inspiral and plunge. 

Let us first recall that a lot of effort has been devoted in recent years to the analytical computing, by means of 
post-Newtonian (PN) expansions in powers of v^/c^ ^ GM/c^r, of the equations of motion, and the GW emission, 
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of comparable-mass binar y s ystems. The equations of motions have been computed to /c^ (3 PN) accuracy by 
two separate groups, |]l^^l5|| and |]l6|-p^, and the two results have been shown to agree ||l^,|2^. Until recently, 
there remained (in both approaches) an ambiguous parameter, ojs, linked to the problem of regularizing some badly 
divergent integrals arising at the 3 PN level. In a recent work, using an improved regularization method (dimensional 
continuation), the first group has succeeded in determining without ambiguity the value of lUs , namely uJs — 
This unique determination of the 3PN equations of motion is consistent with an old argument of 122] showing that it 
should be possible to model black holes by point particles without ambiguity up to the 5 PN level (excluded). The 
emission of GW is unambiguously known to w^/c^ (2.5 PN) accuracy and has recently been formally computed to 

/c^ (3.5 PN) accuracy |2J], modulo the appearance of several ambiguous parameters (^, k, Q) linked to the problem 
of regularizing some divergences arising at the 3PN level. Dimensional regularization is expected to determine without 
ambiguity the values of ^, k, and C, but has not (yet) been applied to the radiation problem. 

We wish to emphasize that such high-order PN results are a necessary^ but not by themselves sufficient, ingredient 
for computing with adequate accuracy the gravitational waveform of coalescing binaries. Indeed, it was emphasized 
long ago ^ that the PN series (written as straightforward Taylor series in powers of some parameter s ~ v/c) 
become slowly convergent in the late stages of binary inspiral. A first attempt was made in [p5| to improve the 
convergence of the PN-expanded equations of motion so as to determine the (crucial) location of the last stable 
(circular) orbit (LSO) for comparable-mass systems. However, further work ^6 27,^ has shown the unreliability (and 
coordinate-dependence) of this attempt. There is, however, no reason of principle preventing the existence of gauge- 
invariant "resummation methods" able to give reliable results near the LSO. Indeed, as emphasized in and |ic| ] 
most coordinate-invariant functions (of some invariant quantity x ^ jt? ^ GM/c^r) that one wishes to consider 
when discussing the dynamics and GW emission of circular orbits are expected to have a singularity only at the 
"light ring" (LR) value of x [last possible unstable circular orbit]. If we trust (for orders of magnitude considerations) 
the small mass-ratio limit {i/ = fi/M = mi 1112/ {mi + 7712)^ 1), we know that xlr — 1/3 is smaller by a factor 
2 than xlso — 1/6. If the functions f{x) we are dealing with are meromorphic functions of x, the location of the 
expected closest singularity (xlr) determines their radius of convergence. Therefore, we expect that, for x < xlr, 
the Taylor expansion of f{x) will converge and will behave essentially like {^/^lrY^- In particular, one expects 
/(a^LSo) ~ X^n (^^LSo/a^LR,)" ~ J2 This heuristic argument suggests a rather slow convergence, but the crucial 

point is to have some convergence, so that the application of suitable resummation methods can be expected to 
accelerate the convergence and to lead to numerically accurate results from the knowledge of only a few terms in 
the Taylor expansion. There exist many types of resummation methods and none of them are of truly universal 
applicability. As a rule, one must know something about the structure of the functions f{x) ~ fo + fi x + f2 x'^ + ■ ■ ■ 
one is trying to resum to be able to devise an efficient resummation method. Refs. [Q, ]^, ]^ and ]|l^ have studied in 
detail the various functions that might be used to discuss the GW fiux and the dynamics of binary systems. This work 
has led to selecting some specific resummation methods, acting on some specific functions. For what concerns the 
GW flux we refer to Fig. 3 of Q] for evidence of the acceleration of convergence (near the LSO) provided by a specific 
resummation method combining a redefinition of the GW flux function with Pade approximants. We wish here, for 
the benefit of the skeptics, to exhibit some of the evidence for the acceleration of convergence (near the LSO) in the 
description of the 2-body dynamics provided, at the 3 PN level, by a resummation method defined by combining |q] 
and ]|l^. Specifically, we mean the combination of the effective-one-body (FOB) approach (further discussed below) 
and of a suitable Pade resummation of the effective radial potential at the n PN level: Ap^ (u) = P^^ [Tn+i{A(u))] (see 
below) . 

Let us consider a sequence of circular orbits, near the LSO and for two non-spinning black holes. In the FOB 
approach the circular orbits are obtained by minimizing a certain effective radial potential, defined by fixing the total 
orbital angular momentum L in the Hamiltonian. The most natural variable defining the one-parameter sequence of 
circular orbits is then simply the angular momentum L. It is therefore natural, for the purpose of this work, to measure 
the separation between the two holes (in a gauge-invariant and approximation-independent) way by conventionally 
defining a £-radius Ri = GM ri , such that the (invariantly defined) total orbital angular momentum L = GiiM£ is 
given by i'^ = rf/{ri — 3), i.e. by the relation holding for a test particle in a Schwarzschild spacetime. [Here, and in 
the following, we shall often set c — 1 and/or G = 1, except in some (final) formulas where it might be illuminating 
to reestablish the dependence on c and/or G.] As the problem is to know whether the resummation method of the 
PN-expanded two-body dynamics is efficient, we compare in Table ^ the total energies frcai of the binary system, 
computed using 1 PN, 2 PN and 3 PN information, for circular orbits at ^-radii ri = 12, 11, 10, 9, 8, 7, and 6. We give 
the values of the binding energy per unit (total) mass, e = {£/M) — 1, for the equal-mass case (mi = 7712; v = 1/4). 

The numbers displayed in Table || illustrate the efficiency of the resummation method advocated in ]^,0 . For 
ri = 12 the fractional difference in binding energy between the 1 PN approximation and the 3 PN one is 0.74%, 
while even for = 6 this difference is only 3.6%. These numbers indicate that, even near the LSO, the Pade- 
improved effective-one-body approach is a rationally sound way of computing the 2-body dynamics. There are no 
signs of numerical unreliability, as there were in the calculations based on the straightforward coordinate-dependent, 
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-1.388 



TABLE I. Binding energies e = {£^cai/M) ~ 1 for circular orbits of equal-mass (non-spinning) binary systems near the 
LSO. The invariant dimensionless ^-radius ri = Ri/GM is defined in the text. The binding energy is computed from the 
Pade-resummed Effective-one-body Hamiltonian at three successive post-Newtonian approximations: 1 PN, 2 PN and 3 PN 
(with LOs =Q )■ 



PN-expanded versions of the equations of motion ||2^ or of the Hamiltonian [g6 27 1 (which gave results differing by 
O (100%) among themselves, and as one changed the PN order). We shall see below that the robustness of the 
PN-resummation exhibited in Table ^ extends to a large domain of the parameter space of spinning black holes. 

As we do not know the exact result, and as current numerical simulations do not give reliable information about the 
late stages of the quasi-circular orbital dynamics of two black holes (see below), the kind of internal consistency check 
exhibited in Table ||is about the only evidence we can set forth at present. [Note that, from a logical point of view, 
the situation here is the same as for numerical simulations: in absence of an exact solution (and of experimental data) 
one can only do internal convergence tests.] Ideally, it would be important to extend the checks of Table | to the 4 PN 
level (to confirm the trend and see a real sign of convergence to a limit) but this seems to be an hopelessly difficult 
task with present analytical means. Finally, let us note that the fact that one can concoct many "bad" ways of using 
the PN-expanded information near the LSO (exhibiting as badly divergent results as wished) is not a valid argument 
against the reliability of the specific resummation technique used in 10 1 and here. An ambiguity problem would 
arise only if one could construct two different resummation methods, both exhibiting an internal "convergence" (as 
the PN order increases) as good as that illustrated in Table |, but yielding very different predictions for physical 
observables near the LSO. This is not the case at present because the comparative study of jl^] (see Table I there) 
has shown that the EOB approach exhibited (when lUs ^ —9) significantly better PN convergence than a panel of 
other invariant resummation methods. 

In the following we take for granted the soundness of the effective-one-body resummation approach and we show how 
to generalize it to the case of two (moderately) spinning black holes. Let us first recall that the basic idea of the EOB 
approach was first developed in the context of the electromagnetically interacting quantum two-body problem [p8|, 
pof (see also ||30[l ) . A first attempt to deal with the gravitationally interacting two-body problem (at the 1 PN leve ) 
was made in |31| (see also ^^ ). A renewed EOB approach (which significantly differs from the general framework set 
up by Todorov and coworkers |2^j3l[]) was introduced in The latter reference showed how to apply this method at 
the 2 PN level. It was then used to study the transition between the inspiral and the plunge for comparable masses, 
and, in particular, to construct a complete waveform covering the inspiral, the plunge and the final merger (see ]3^ ] 
for the physical consequences of this waveform) . More recently, Ref. showed how to extend the EOB approach to 
the 3 PN level (this required a non trivial generalization of the basic idea). 

Before entering the details of our way of introducing the spin degrees of freedom in the EOB approach, let us state 
our general view of the usefulness of the EOB method in this context. As we shall discuss below, the present work 
(which incorporates spin effects at leading PN order) can only be expected to give physically reliable results in the case 
of moderate spins ( flp ^ 0.3, see below). However, the EOB approach, far from being a rigid structure, is extremely 
ficxible. One can modify the basic functions (such as A{u)) determining the EOB dynamics by introducing new 
parameters corresponding to (yet) uncalculated higher PN effects. [ These terms become important only for orbits 
closer than 6GM, and/or for fast-spinning holes.] Therefore, when either higher-accuracy analytical calculations 
are performed or numerical relativity becomes able to give physically relevant data about the interaction of (fast- 
spinning) black holes, we expect that it will be possible to complete the current EOB Hamiltonian so as to incorporate 
this information. As the parameter space of two spinning black holes (with arbitrarily oriented spins) is very large, 
numerical relativity will never be able, by itself, to cover it densely. We think, however, that a suitable "numerically 
fitted" (and, if possible, "analytically extended") EOB Hamiltonian should be able to fit the needs of upcoming GW 
detectors. The present work should be viewed as a first step in this direction. 

The present paper is organized as follows. In Section || we show how to incorporate (in some approximation) the 
spin degrees of freedom of each black hole within a 3 PN-level, resummed effective one-body approach. In Section III 
we study some of the predictions of our resummed dynamics, notably for what concerns the location of the transition 
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between the inspiral and the plunge. Section IV contains our concluding remarks 



II. EFFECTIVE ONE-BODY APPROACH, EFFECTIVE SPIN AND A DEFORMED KERR METRIC 

A. Effective one-body approach 

Let us recall the basic set up of the effective one-body (EOB) approach. One starts from the (PN-expanded) two- 
body equations of motion, which depend on the dynamical variables of two particles. One separates the equations 
of motion in a "conservative part" , and a "radiation reaction part" . Though this separation is not well-defined at 
the exact (general relativistic) level it is not ambiguous at the 3 PN level (in the conservative part) which we shall 
consider here |^ We shall henceforth consider only the conservative part of the dynamics. [We leave to future work 
the_^eneralization to spinning black holes of the definition of resummed radiation reaction effects which was achieved 
for non-spinning black holes.] It has been explicitly shown that the 3 PN^ dynamics is Poincare invariant [ p^ , 
The ten first integrals associated to the ten generators of the Poincare group were constructed in |jl^ (sec also 
|20|). In particular, we have the "center of mass" vectorial constant K = G — t P. This constant allows one to define 
the center of mass frame, in which K = 0, which implies P = and G = 0. We can then reduce the PN-expanded 
two-body dynamics to a PN-expanded one-body dynamics by considering the relative motion in the center of mass 
frame. This reduction leads to a great simplification of the dynamics. 

Indeed, the full 3 PN Hamiltonian in an arbitrary reference frame [|l5| contains O(IOO) terms, while its center-of- 
mass-reduced version contains only 24 terms. However, this simplification is, by itself, insufficient for helping in any 
way the crucial problem of the slow convergence of the PN expansion. One should also mention that the use of an 



Hamiltonian framework (like the ADM formalism used in |12| , |l3| , [l0 15 2l|]) is extremely convenient (much more so 
than an approach based on the harmonic-coordinates equations of motion, as in [^6]-^ 20 1). Indeed, on the one hand 
it simplifies very much the reduction to the center-of-mass relative dynamics (which is trivially obtained by settin; 
p^gj = Pi — —P2)i s-iid on the other hand it yields directly (without guesswork) an action principle for the dynamic; 
We shall find also below that an Hamiltonian approach is very convenient for dealing with the spin degrees of freedom. 

Up to now, we only mentionned the dynamics of the orbital degrees of freedom, i.e. (in the order-reduced Hamilto- 
nian formalism) the (ADM-coordinate) positions and momenta Xi,X2,Pi, P2 of the two black holes^. After reduction 
to the center-of-mass frame {P — Pi + P2 = 0), and to the relative dynamics {x = Xi — X2, p = Pi = —P2), one ends 
up with a canonical pair x, p of phase-space variables. 

The addition of spin degrees of freedom on each black hole is, a priori, a rather complicated matter. If one wished 
to have a relativistically covariant description of the dynamics of two spinning objects, one would need not only to 
add, in Einstein's equations, extra (covariant) source terms proportional to suitable derivatives of delta-functions 
(spin dipoles), but also to enlarge the two-body action principle to incorporate the spin variables. The first task is 
doable, and has been performed (to the lowest order) in several works, such as Refs. ]3^-|3^]. For other works on the 
relativistic equations of motion of black holes or extended bodies (endowed with spin and higher multipole moments) 
see ]39|-^. Recently, Ref. has tackled the next-to-leading order contribution to spin-orbit effects. [We consider 
here only the case of interacting, comparable massive objects. The problem of a spinning test particle in an external 
field is simpler and has been dealt with by many authors, such as Mathisson |4^, Papapetrou |^^, etc. . .] On the 
other hand, the second task (incorporating the spin degrees of freedom in an action principle) is quite intricate. First, 
it has been found that, within a relativistically covariant set up for a spinning particle, the Lagrangian describing the 
orbital motion could not, even at lowest order in the spin, be taken as an ordinary Lagrangian L(x, i), but needed to 
be a higher-order ) Wm. Second, a relativistically covariant treatment of the spin degrees of freedom 



We expect real ambiguities to arise only at the /c^° ~ 5 PN level, because this corresponds to the square of the leading, 
/c" ~ 2.5 PN, radiation reaction terms. 
^Henceforth, '3 FN' will mean the conservative 3 PN dynamics, i.e. 7V + 1PN-I-2PN-I-3 PN. 

^Note that a subtlety arises at 3 PN in that the Hamilton action principle involves derivatives of the phase space variables. 
However, it was shown in jl^] how to reduce the problem to an ordinary Hamiltonian dynamics by means of a suitable 0{v^/c^) 
shift of phase-space variables. We henceforth assume that we work with the shifted variables defined in [ p^ . 

*We recall that the high-order perturbative, PN-expanded, calculations of the dynamics of two non-spinning compact objects 
model these objects by delta-function (monopole) sources. The supports of these delta-functions define the coordinate "posi- 
tions" of the compact objects. As explained in [E2| these "positions" physically correspond to some "centers of the gravitational 
field" generated by the objects. 
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is an intricate matter, involving all the subtleties of constrained dynamical systems, even in the simplest case of a 
free relativistic top |Q. Contrary to the case of the spin- independent EOB where it was easy to use the test- mass 
results to constrain the EOB Hamiltonian, the Mathisson-Papapetrou dynamics of spinning test masses in external 
gravitational fields is rather complicated and cannot easily be used to constrain the spin-dependent EOB Hamiltonian. 
(It might, however, be interesting to try to do so.) Fortunately, there is a technically much lighter approach which 
bypasses these problems and simplifies both the description of orbital degrees of freedom and that of spin degrees of 
freedom. This approach is not manifestly relativistically covariant. This lack of manifest Poincare covariance is not 
(in principle) a problem at all for two reasons: (i) it does not prevent the expected global Poincare covariance of 
the two-body dynamics to be realized as a phase-space symmetry (as was explicitly proven, at 3 PN, for the orbital 
degrees of motion in Ref. p^]), and (ii) as we are, at this stage, mainly interested in the description of the relative 
motion in a specific (center-of-mass) Lorentz frame, there is no physical need to enforce any boost invariance. This 
non covariant approach to the gravitational interaction of spinning objects stems from the classic work of Breit on 
the electromagnetic interaction of (quantum) spinning electrons (see, e.g. p6|| ) and has been developed in a sequence 
of papers |47'50|. [The change of variables needed to pass from the covariant, higher-order Lagrangian description to 
the non-covariant, ordinary Lagrangian has been discussed in several papers, e.g. ||45| , |5l| , ^ .] In the present paper, we 
shall combine the non manifestely covariant, ADM-Hamiltonian treatment of the orbital degrees of freedom of ||l2|,|l5|, 
with the similarly non covariant, but Hamiltonian, treatment of the spin degrees of freedom of 49 1. Moreover, as 
is explained below, we shall improve upon ||4^-|50t] in using a direct Poisson-bracket treatment of the dynamical spin 
variables. 

Finally, our starting point (for the effective one-body approach) is a PN- expanded Hamiltonian for the relative 
motion of two spinning objects of the form 



H^Zix,p,Si,S2) 



H^^{x,p, Si, 82) + H^sssi^^P^ Si,S2) 



(2.1) 



Here, -ff™ denotes the PN-expanded orbital Hamiltonian, which is the sum of the free Hamiltonian Hq 



pI + \frn\ 



+ P2 ^"^^ the monopolar interaction Hamiltonian generated by the source terms 



proportional to the masses. 



iJ;™^ ~ mi 7712 +'Tii 7712 -fmi m| -t-m^ m^+rri^ rr3^^m\ m2 + 
It is explicitly known up to the 3 PN level (i.e 



Before the reduction to the center-of-mass frame i?;™^ has the symbolic structure: 

nil m2{l + mi+m2 + {mi+m2)^ + {mi+m2)^ + - ■ •). 
up to velocity-independent terms oc mi 7712(7711 -|- m2)'^). Af- 



ter reduction to the center-of-mass frame the PN expansion of H^^^ reads (with M = mi -|- m2, /U = mim2/M, 
V = [ijM = mi m2/(mi 4- m2)^, r = \x\,p = pl [i, x = x/GM) 

Hl^{x,p) = Mc^ + Hn{x,p)+Hipj^{x,p)+ H2pn{x,p) + H3Pj^{x,p) , (2.2) 



Hn{x,p) 



El 



1^2 1 
2^ 



(2.3) 



Hipf^{x,p) = — 



I (3j/ - 1) p** - i [(3 + i^)p^ + v{n ■ p) 



1 
2F 



(2.4) 



2PN 



{x,p) - ^ 



bv + bv'^)p + 



1 , ^ 1 ' 

7(1 + 3-)-^ 



(2.5) 



^?3PN(a;,p) 



128 



(-5-l-35i/-70i/2-f 35i/3)p** 



109 
"T2" 



21 
32 





1 






f4 



(2.6) 



We have exhibited (for illustration) in Eqs. ( |2.5|) and (^^) only the first and the last terms. We refer to jS^] 
for the full (center-of-mass) 2 PN Hamiltonian (7 terms in all), and to jl^Jl^,|l^ for the full (center-of-mass) 3 PN 
Hamiltonian (11 terms in all). Our effective one-body treatment will take into account the full 2 PN and 3 PN 
structures, but in a very streamlined way which will be explicitly displayed below. 
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The other terms in Eq. (2.1) denote the various spin-dependent contributions to the Hamiltonian: respectively the 
terms hnear (_ff™), quadratic (^55), cubic (Hggg), etc. . . in the spins ^i, ^2. Before reduction to the center-of-mass 
frame they have the symboUc structure: 



H 



PN 



Si TO2(1 + mi + 7712 + {mi + 7772)^ + • • •) 

+ 'S'2 777l(l + 7771 + 777 2 + ("ll + '772)^ + ' ' ') 



(2.7) 



H 



Is ~ S"! 7772 I — + 1 + mi + 7712 + •••)+ 5*1 5*2(1 + 77ll + 7772 + • • •) 



1 

7771 

1 
7772 



5| 7771 ( — h 1 + 7711 + 7772 



(2.8) 



etc. . . We shall explain below the occurence of the terms quadratic in the spins and inversely proportional to a mass. 
In contradistinction with the case of the orbital Hamiltonian which has been worked out with a high PN accuracy, 
only the simplest spin-dependent terms have been explicitly derived, namely the lowest PN-order term in -ff™, whose 
center-of-mass reduction reads M 



H^''{x,p,Si,S2) 



2G 



S m2 , 
1^1^^' 



4 7712 



(2.9) 



and the lowest PN-order one-graviton exchange contribution t o t he bilinear term (cx 5*1 S'2) in H^g . Ref. |42| ] 
contains some information about the C'(^) corrections in Eq. ( ^.9[ ), but, because of the use of different gauge and 
spin conditions, not in a form which can be directly used to derive these corrections. We shall discuss the spin-bilinear 
contribution (oc Si S2) below, together with the leading spin-quadratic contributions oc Si 7772/7771 -I- 5*1 7711/7712. 

Before going further, let us make clear that, before and after any type of resummation, the dynamics entailed by 
the Hamiltonians we shall consider H{x,p, Si, S2) follow, for all degrees of freedom, from the basic Poisson brackets 



(2.10) 



{SI Si} = e'^'^S'^., 



{S^2,Si} = e'^'S^. 



(2.11) 
(2.12) 



= {x\x^} = = {SI, Si} - {x\Si} = K 



• Si} 



{pi,si} = {pi,si}. 



The (real) time evolution of any dynamical quantity f(x,p. Si, S2) is given by 



(2.13) 



dt 



f{x,p, 81,82) ^ {/,i?roal}, 



(2.14) 



where the Poisson bracket (PB) {/, i?rcai} is computed from the basic PB's (|]i^-(|t|) by using the standard PB 
properties (skew symmetry: {/,<?} = — {ff, /}i Leibniz rule: {f,gh} = {f, q} h + g {f, h }, an d the J acobi identity: 
{f,{g,h}} + {g,{h,f}} + {h,{f,g}} — 0). The simplest way to prove the statements ( ^.1(]| )-( 2.14 ) is to consider 
our dynamics as the classical limit of the quantum dynamics of a system of gravitationally interacting spinning 
particles. Surprisingly, though Refs. |47|4^ derived (a la Breit) the spin-dependent contributions to the Hamiltonian 
by a quantum ro ute, they never noticed that they could very simply derive the spin evolution equations by using 
the PB's ( ^.llD , (b.l2D . They had to go back to a Lagrangian formalism and add some explicit spin kinetic energy 



terms 



/2i£i) 



to derive the spin evolution equations. Note also that we have kept the label "real" on the 

Hamiltonian in Eq. ( |2.14| ) to distinguish the evolution with respect to the real time (associated with the original two- 
body system) from the evolution generated by the effective Hamiltonian to be introduced below (which is associated 
with an auxiliary, effective time). 

Before generalizing it, by including the spin degrees of freedom, let us recall the results of j|] (2 PN level) and 
I pot (3 PN level) concerning the effective one-body "upgrading" of the PN-expanded orbital Hamiltonian iJ™(a;,p). 
Again, the simplest way to motivate it is to think of our dynamics as the classical limit of a quantum dynamics 
defined by some hermitian Hamiltonian operator _fforb(3', — iftV). We are mainly interested in the bound states of 



6 



i/orb- It is crucial to note that the orbital Hamiltonian (2.2)-p^) is symmetric under an 0(3) group (corresponding 



to arbitrary rotations of the relative position x — Xi — X2, in the center-of-mass frame). Therefore the quantum (and 
classical) bound states will be labelled (besides parity) by only two quantum numbers: (i) the total orbital angular 
momentum L'^ = L(L + h), and (ii) some "principal quantum number" N (such that {N — L)/h counts the number 
of nodes in the radial relative wave function). Both L and N are quantized in units of h. The full list of two-body 
bound states is thereby encoded in the formula giving the bound state energy as a function of the two quantum 
numbers L and : E.^^i = E,^^i{L,N) = Mc^ - \ ^i{G rrn 1112? / N"^ + Ei^L, N) + E2p^{L, N) + E3pt^{L,N). 
The basic idea of the effective one-body method is to map (in a one-to-one manner) the discrete set of real two-body 
bound states Ereed{L, N) onto the discrete set of bound states of an auxiliary ("effective") one-body Hamiltonian 
cff,Pcff)- Because of the special labelling by (only) two integer quantum numbers L/h^ N/h one is naturally 
led to imposing that: (i) the EOB Hamiltonian be spherically symmetric^, and (ii) the integer valued quantum 
numbers be identified in the two problems, i.e. L/h — Lcff/h and N/h = Neff/h. On the other hand, one can 
(and one a priori should) leave free a (one-to-one) continuous function / mapping the real energies onto the effective 
ones: Ecs{L,N) = f{Ei-cii\{L,N)). Evidently, for this method to buy us anything we wish the effective dynamics 
to be significantly simpler than the original Hrca.i{x,p), and, in particular, to reduce, in some approximation, to the 
paradigm of the simplest gravitational one-body problem, namely the dynamics of an (effective) test particle in some 
(to be determined) effective metric g'r^{x^ff). Remarkably enou gh, it was found in that such a mapping between 



the very complicated real two-body orbital 2 PN Hamiltonian (|2.2| )-(2.5) and the usual ("geodesic") dynamics of a 



test particle of mass fi = mi m2/(mi -I- 7712) in a very simple (spherically symmetric) effective metric 

dsls = -A(r,ff) dtl^ + ^^p^ drl^ + rlside^ + sin^ 9 d^') , (2.15) 



is possible, if and only if the energy mapping E^s = /(-Ercai) is given by 



fxc^ 2 mi TO2 



(2.16) 



Remarkably, the simple energy map (2.16) (which is here determined by our requirements) coincides with the energy 
map introduced in several other investigations |2^ , 0| (and is simply related to the one defined a . pri ori in ||2^,^,| 



Recently, the problem of mapping the extremely complicated real two-body 3 PN Hamiltonian (2.6) onto an effective 
one-body dynamics has been solved jl^. Again the result is remarkably simple, though less simple than at the 2 PN 
level. Indeed, it was fount^ that the effective one-body dynamics was given by an Hamilton- Jacobi equation of the 
form (with — dS/d 



cSJ 



= f^'+ 9:^pfpf + Q4(p'"^) , (2.17) 



with a simple (spherically symmetric) effective metric of the form ( 2.15 ) and some additional quartic-in- momenta 
contribution Q/^lp). Remarkably, it was found that, at the 3 PN level, the energy mapping is again uniquely determined 
to be the simple relation ( 2.16 ). As for the metric coefficients of the (covariant) effective metric 5°^, and the quartic 



terms Q^lp), they were foimd to be 

A{r) ^l-2u + 2i^u^ + a4{i^)u^ , (2.18) 

D{r) = 1~ 6iyu'^ + 2(3v -26) lyu^, (2.19) 

9^^^^2{A-3iy)iyu^n-p)\ (2.20) 



^We are making this very explicit because some people, when they hear about the EOB approach, think that the effective 
metric describing the one-body dynamics should, at some level of approximation, include some Kerr-like features to model the 
velocity-dependent two-body interactions. This is not true for the orbital dynamics, whatever be the PN accuracy level. On 
the other hand, we shall see that we need Kerr-like features to accomodate the intrinsic spin effects. 

®In fact, found that it was possible to map the real dynamics onto the geodesic dynamics of a test particle. However, 
both the effective metric and the modified energy map needed for this representation are rather complicated. It was felt that it 
is more convincing to keep a simple effective metric, and a simple energy map, but to relax the constraint of geodesic motion. 
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where u = GM/r, P ^ p/lJ- and 
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41 
32 



(2.21) 



As we said above, the correct value of Ug has been recently found to be simply ojs = 0. We shall find, however, 
convenient to keep lOs as a free parameter in order to assess the quantitative importance of 3PN effects. 

Let us emphasize again the streamlined nature of the effective one-body description of the orbital dynamics. Succes- 
sively, as the PN order increases, one can say that: (i) the 6 terms of the Newtonian plus first post-Newtonian relative 
Hamiltonian (2^), (2^) can be mapped (via ( ^.16 ) and a canonical transformation of {x,p)) onto geodesic motion in 
a Schwarzschild spacctimc of mass M (i.e. pn = 1 — 2 u, Di pn = 1), (ii) to take into account the 7 additional terms 
entering -ff2PN, Eq. (|2.5|), it is enough to add +2iyifl to A{r) and — 6 z/u^ to D{r), and (iii) to take into account the 



11 additional terms entering j?3 PNi it is enough to further add -I- 4(1^) u to A{r) and -I- 2(31^ — 26) to D{r), and 



to add the simple quartic term ( 2.2C ) to the mass-shell condition ( 2.17 ). Note that the effective one-body dynamics 
is a "deformation" of a geodesic dynamics for a particle of mass /i in a Schwarzschild spacetime of mass M , with 
deformation parameter the symmetric mass ratio i' — n/M — mim2/(mi -I- 7712)^ (0 < z^ < ^, with 1/ being small 

have not yet introduced any particular resummation technique. The effective 



if 77ii <^ 7772 or 7772 <C 777i, and reaching its maximal value of i when ttii — 1712). Note also t h at, at this stage, we 



quantities ( 2.18|) -( ^.20 ) are still given 



as straightforward PN expansions in p o wers of w = GM/c r. However, this already means an appreciable gain over 



the original PN expansions, Eqs. (2^)-(2^). Indeed, there are far less terms in the "effective" PN expansions and 
they generically have significantly smaller coefficients (which are now all multiplied by v < 1/4). For instance, the 
"radial" potential determining the circular orbits is now fully encoded in the simple function (^qq = —A[rcg) which 
differs from the test-mass (Schwarzschild) result, ^(^^0) — 1 — 2 72, only by two numerically smallish terms when 
one is above the last stable orbit. Indeed, when > QGM/c^, 2vu^ < 0.23% and 04(7^) < 0.36% (for ujs = 0). 
When working at the 2 PN level it was, in fact, found unnecessary [|j to further "resum" the effective PN expansions 
of A{r) and D{r). As they stand, they led to small deviations from the test-mass {ly 0) results. However, it was 
found in fl^ that the situation is not quite as rosy at the 3 PN le vel. Because of the largish coefficient in Eq. ( p. 21 ), 
T " M^^^ 18.688, the additional term -1-04(7^)72* in Eq. ( 2.18 ) significantly modifies the qualitative behaviour of 
the metric coefficient A{r) for < 6 GM. In particular, when LOg — 0, the straightforward PN-expanded function 
A{r) no longer features a zero near r^s = 2 GM for all possible values of the deformation parameter < < i. As 
this zero (which corresponds to the Schwarzschild horizon) is a crucial qualitative feature of the — limit, it was 
argued in ||l^ that one should Pade-resum the PN-expansion of A{r) so as to ensure the stable presence of a similar 
"horizon" when < i' < j- We shall also do so here, but only after having introduced the spin effects, which modify 
the radial function which is the analog of A{r). 



B. Effective spin 



Let us now tackle the central task of this work: to introduce spin effects in the effective one-body approach. Let us 
first emphasize that the ambition of the present work is somewhat limited. Our main goal is to derive a spin-dependent 
EOB Hamiltonian which is physically reliable for small and moderate spins. We shall proceed toward this goal in 
successive steps. First, we consider only (at the lowest PN order where they enter) the interactions which are linear 
in the spins. Then, we shall incorporate (to some approximation) the interactions which are quadratic in the spins. 
Contrary to the case of the spin-independent interactions where many years of work have yielded high-PN-accuracy 
information which allowed one to refine and test the EOB approach, one does not have in hand enough information 
for gauging the reliability of the spin-dependent interactions in the case of fast spins. As a consequence, we will be 
able to trust the presently introduced EOB Hamiltonian only when spin effects are not too large. As we said in the 
Introduction, one will need new information (either from numerical relativity, or from improved analytical methods) 
to find a reliable form of the EO B H amiltonian for large spins. 

We have written down in Eq. (2.9) above the contribution to the real, PN-expanded, two-body Hamiltonian which 
is linear in the two spins. Our first proposal is to map this contribution to the spin-orbit coupling of a (spinless) 
effective particle moving in a suitably "spinning" effective metric, i.e. some type of generalized Kerr metric. If we 
formally consider Si and S2 as deformation parameters (on top of the basic orbital deformation parameter i/), the 
effective dynamics we are looking for should be a "spin deformation" of the currently most accurate orbital dynamics, 
as described by the 3 PN effective dynamics Eqs. ( 2.15D -(2.21). In particular, we should keep the non-geodesic 
Hamilton- Jacobi equation ( 2.17 ). Note that the spin effects (notably Eq. (^^)) break the 0(3) symmetry of the 
orbital interaction. At the quantum level, this means that spin interactions lift the degeneracy {— lack of dependence 
on the "magnetic" quantum number Lz/h) of the orbital energy states. This shows that the effective (co)metric g'^^ 
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entering Eq. (2.17) should no longer be spherically symmetric, but should contain special directions linked to the 

(2.22) 



spins. For generality, let us first consider an arbitrary (time-independent) effective cometric 

9"s Pa Pf3 = ffeff + 2 9oS Po Pi + gls Pi Pj ■ 



Let us define 



I.e. 



00 ' ' — 

-'off 



Oi Oj 
t/cff ffcff 



9oS 



2 ' 9cS 



P ij ij P P'^ 

-— , 5off = 7^ — 



(2.23) 



(2.24) 



The effective energy i?off = ~P°q^ is conserved (because of the assumed stationarity of g^jf ). Using the parametriza- 
tion (|2.23|), Eq. (|2.17|) reads 



(2.25) 



Solving Eq. ( ^.25 ) for Ecs (using the fact that Qi{p), Eq. ( 2.20| ), depends only on the pis) yields the effective 
Hamiltonian 



(2.26) 



where we have suppressed, for readability, the labels "eff" on the orbital phase space variables a;*^^, The ellipsis 
in the arguments of H^s are added to remind us that H^s will ultimately also depend on the spin variables Si, S2, 
which enter the metric coefficients a, 7'-' as parameters. 

We also assume that Eq. (2.16) (which was found to hold at 1 PN, 2 PN and 3 PN) still holds. Solving it for the 
real energy E'rcai in terms of the effective one finally yields the real Hamiltonian 



Ere^l = i?rcal(a5, P, . . .) = M J 1 + 2 U 



HcS — P- 



(2.27) 



We recall that, at the linearized level and at the lowest PN order, the addition of a spin S'off onto an initially 
spherical symmetric metric leads to an off-diagonal term in the metric: 



(2.28) 



Inserting this term in ( p. 26 ), and expanding (2.27) in a PN series yields, as leading spin-orbit coupling (linearized 
in S'off and taken to formal order 0{\/c^)) in -ffrcai, the term 



2G 

5s,it -ffrcal ^ Pt ^ e^^^ p, S^ff x'' . 

This term can exactly reproduce the leading^ two-body spin-orbit coupling (|]^) if we define 



-,off 



(Ti Si + a2 S2 , 



with 



3 1712 3 mi 

cri = l + - , (72 = 1 + - 

4 mi 4 1712 



(2.29) 



(2.30) 



(2.31) 



'^We use here the fact that the real phase-space coordinates a;"^"^', differ only by 0{l/c^) from the effective ones entering 
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C. A deformed Kerr metric 



Remembering that the main message of the effective one-body method is that the orbital dynamics of two 
comparable-mass black holes can be described in terms of a slightly deformed (with deformation parameter i') 
Schwarzschild metric, we expect that the orbital-plus-spin dynamics of two black holes can be described in terms 
of some deformation of the Kerr metric. In other words, we are expecting that not only the effects linear in the 



spins, such as Eq. (2.9), but also the spin-dependent non-linear effects, can be described in terms of some deformed, 
effective Kerr metric. At this stage it is th erefor e very natural to construct a suitable "deformed Kerr metric" which 
combines the orbital deformations (2. IS), (2.19) with the full spin effects linked to the "effective spin" (2.30). After 
constructing this deformed Kerr metric, we shall a posteriori check that it approximately incorporates the expected 
two-body interactions which are quadratic in the spins. 

Let us start from the simplest form of the Kerr cometric, underlying its separability properties |53| 



q/3 



+ cos^ 9 



\— (p^ + asin^ 9pt)'^ 



AK{r) 



{{r^ +a'^)pt+ap^)' 



(2.32) 



with A/f (r) — r"^ — 2Mr + a^. In the non-spinning limit (a 0) the coefficients of an d p^ become, respectively, 
A/f (r')/r^ and — r^/A/f (r). However, we know that in this limit we should get (f rom ( 2.15 )) A{r)/D{r) and —1/A{r), 
respectively. It is therefore very natural to generalize the Kerr metric ( |2.32 ) (while still keepi ng its separability 
properties) by assuming that the coefficients of the first and last terms in the square brackets of ( |2.32| ) involve two 
different functions of r, say A, (r) and — 1/A((r), whose product reduces to —1/D{r) when a — > 0. This reasoning 
leads us, as sim plest^ possibility for combining spin effects with orbital effects, to postulating that the effective metric 
entering ( 2.17| ) has the form 



a/3 



1 



r'^ + cos^ 9 
Ar{r)pl+pl + 



sm 



<Pt) 



Mr) 



a^)Pt 



(2.33) 



with 



AP(r) A^^"(r) + 



(2.34) 



Here, the superscripts "PN" indicate that, at this stage, we are only comparing PN expansions. We already know from 
the 3 PN study of |l^] that this is unsatisfactory because it tends to change the qualitative behaviour of t he ra dial 
functions, and, in particular, the presence of a horizon in the metric (2.33). To get a regular horizon in Eq. (2.33) we 
need the two functions At(r) and Ar(r) to have a zero at the same value of r. The simplest (and most robust) way 
of ensuring this is (as discussed in [fol) to define them as 



Mr) 



r'Pl 



A'''\r) + 



A,(r) = At{r) 



D{t) 



PN 



(2.35) 



Here, with u = l/r, denotes the Nn{u)/Dm{u) Pade of a certain PN-expanded function f (u 



Co + Ciu + C2U" + ■ ■ ■ + Cn+m m"+™ {Nn{u) and Dm{u) being polynomials in u of degrees n and m, respectively). We 
do not write down the (uniquely defined) explicit expression of 



A{u) = P3 [A^^^ (m) + u^] = P3 [1 - 2u + + 2v^i:^ + ai{v) 

1 + ni u 



1 



diu -\- d2'u?' 



ds 



We leave untouched the dependence on a to ensure that, when GM with a being fixed, the metric g^g be Minkowski in 
disguise. 
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(where u = GM/r, a = a/GM) because: (i) it is rather compUcated and not iUuminating, and (ii) modern algebraic 
manipulators compute it directly from its Pade definition. 

In the definition of Ar(r) (which is less important than that of Af(r)) we have factorized the Padeed Af(r) and 
assumed that it was enough to work with the non-rcsummcd PN-expansion of the inverse of the _D-function, i.e. (from 



( P-19D ) 

{D'^{r)f^ EEl + Qvu"^ + 2{2&~iv)vu^ . (2.36) 

If the need arises, it would be easy to define improved (resummed) versions of D^^{r). Because of the positive 
coefficients in (2.36) the present definition does not interfere (as would the consideration of {D^^{ry)~^) with the 



desired feature of having a simple zero in (r) located at the same value as that in Aj (r) . 

Finally, we shall see later that there are some advantages in defining the quartic- in- momenta contribution Qi{p) in 
the following (deformed) way 

Q4b = — 2 4-3Z. / 2fl ^-P ■ 2.37 

/u^ + cos^ 



Eq. (2.33) defin es q"J ^ only with respect to some (instantaneous) polar coordinate system with 2;-axis aligned with 



the effective spin ( 2.3C ). Such a coordinate system cannot be used for describing the evolution of two gravitationally 
interacting spinning black holes. Indeed, we expect (and shall check below) that the total real Hamiltonian imposes 
some type of precession motion for Si, S2 and therefore ScS- To get the full dynamics of the system we need to rewrite 
g"g in a general, cartesian-like coordinate system. This is achieved by explicitly introducing, besides = x'^/r, the 
quantities {S,s = {Sat S^^ 3'^^^)^/^) 

- " - "^^^^ cos = n's\ p^ = + cos^ 9 . (2.38) 



This leads to the following, cartesian-like, effective metric 



a + 0^ ~ At 



(s X xY , 



(2.40) 



p2 g^J^ = + {5'^ - n^ ) - (s x x)' {s x xY , (2.41) 



At 

from which follows 



P"^' (2.42) 



(r2 -I- a2)2 _ a? At sin^ 9 ' 



(3^ 



a [r2 + a2 - At] (s x xf 



(r2 



- a? At sin2 



(2.43) 



(3' 



(2.44) 



Note that near the "horizon", i.e. as A^ 0, the quantity a ten ds to zero like \/Ai, while and 7*-' have finite 
limits. [The singular last term on the right-hand side of Eq. (2.41) is cancelled near the horizon by the contribution 
+P^P^/a^ to Eq. ^M).] 



Finally, the spin-dependent, real two-body Hamiltonian H^ca\{x, p, Si, S2) is defined by (pi = Pi/ fJ-, Up = GM/ a 
ri* = x'^/r) 
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where we recall that the basic effective Kerr spin vector is defined by 

Mas' EE Slfi EE (Jl S\+(J2Sl., 



(2.45) 



(2.46) 




3l|). The phase-space coordinates appearing in this Hamiltonian are the effective 
by 0(l/c^) terms from the coordinates used in usual PN calculations, such as 



with cTi and f72 defined in E< 

ones (a;*ff ,Pj*^). They differ 

ADM o nes. The evolution equations defined by the Hamiltonian (2.45) ar e obt ained by the Poisson bracket equations 
( ^.10 )~( 2.13| ). Before discussing them let us show how the Hamiltonian ( 2.45| ) contains spin-quadratic effects of the 
good sign and magnitude. 



D. Effects quadratic in the spins 

Note first that if we introduce the "non relativistic" effective Hamiltonian = Hcs — ^c^, and similarly 

ifNREEi7„ai-Afc2, one has 

if,T = « (l + I ff$) . (2.47) 

Therefore, if we are interested in the leading PN approximation to any additional term in HfJ^^, one can neglect the 
(0(l/c^) smaller) difference between iJoff and i?rcai- By this argument, the leading PN approximation to the term 
linear in the spins is 

(r^ + a-^Y — a'^iAt sm 

We write it explicitly in the form in which it appears in our Hamiltonian for the reader to see how the term ( ^.29| ) is 
generated. [The important feature here being that -I- — At ~ 2 GM r at the leading PN approximation.] 

Let us now consider the interact ion te rms in i/rcai or _ffoff which are quadratic in the spins, and therefore quadratic 



in the Kerr- like parameter a, Eq . ( 2.46 ). First, one should remember that most terms of order a^, as they appear in 
the effective metric (2.39)-(2.41), do not directly correspond to physical effects proportional to S^q. Indeed, we are 
using here Boyer-Lindquist-type coordinates which differ, even in the flat space limit GAI — > 0, from usual (flat-space, 
cartesian-like) coordinates by terms of order O(a^) . As we are interested in the leading PN effects quadratic in the 



spins, we can view the Kerr-like metric (2.39)-( ^.41 ) as a deformation, by the a-dependent terms, of the Schwarzschild 



metric (which is the leading PN version of the orbital effective metric). We then expect that the leading physical 
effects quadratic a will be those linked to the a-dependent quadrupole moment deformation of the Schwarzschild 
metric. The quadrupole moment of the Kerr metric (which coincides with our metric when we neglect additional 
2 PN fractional corrections) has been determined to be 

= -M Sj Sj + - Ma^ Sij . (2.49) 
3 

This corresponds to an additional term in the interaction Hamihonian equal to 

^H^Q^^c^-^fi Q„ a„ J - +^ M Ma^ % ^ . (2.50) 
In terms of the effective spin this reads (in standard units) 

HT'^+^o^T;fSlsSi^d,,l. (2.51) 



Such is the prediction from our Hamiltonian. Let us now compare it to the expected real two-body, spin-quadratic 
effects. As sketched in Eq. (Oh there are several sources of spin-quadratic effects. At leading PN order, it is enough to 
consider: (i) the term cx 7712 Sj/toi which arises because of the interaction of the monopole m2 with the spin-induced 
quadrupole moment of the spinning black hole of mass mi, (ii) the term oc nii 5*1/7712 obtained by exchanging 1 — > 2, 
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and (iii) the term cx 5*1 5*2 coming from t he dir ect, one-graviton interaction between the two spin-dipoles. The first 
term is obtained by relabeUing the resuh ( 2.51 ) by /i — > ni2, M — > rrii, Slg — > SI- Therefore the sum of (i) and (ii) 
reads 



HsiSi + HS2S2 



2 \mi 



1 



The term (iii) has been computed in and reads 



Hs,s. ^ SI Si (3 - 5^= ) = S\ Si 8,, - 



It is easily checked that the sum of (2.52) and (2.53), say Hss = Hs^Si + HS2S2 + can be written as 



lGjj_ 
2^Tl 



Hss — n ~ 'VT SqSq d^j - , 



(2.52) 



(2.53) 



(2.54) 



with Sq/M = = a\+ = S\/mi + S^jmi-, i.e. exphcitly. 



1 



m2 



(2.55) 



The resuh (2.54), (2.55) is remarkably similar to the prediction ( ^.5l[ ) (with (2.30), (2.31)). The only discrepancy is 
a 25% difference in the coefficient of the mass ratios in the definition of the effective spin. Though there might be 
physical situations where this smallish difference might play a significant role, we think that in most cases where one 
will be entitled to trust the approximate spin-dependent EOB Hamiltonian introduced here this difference will not 
matter. Indeed, because of the partially ad hoc way in which we constructed our deformed Kerr metric, we cannot 
trust our predictions beyond the domain where spin effects are moderate corrections to orbital effects. However, it 
is useful to incorporate in a qualitatively correct manner the non- linear spin e ffects . This is w hat our prescription 
achieves. For instance: (i) in the limit where, say, mi ^ mi (and |5'2| < mg ) (2.51) and (2.54) become equivalent, 
or (ii) in the case where S^i and ^2 are parallel (in the same direction), (2.51) and (2.54) differ only by a numerical 
factor which is near one for all mass ratios. 

It might, however, be useful to define another Hamiltonian, say H^.^^^, which: (a) reduces (like TJrcai) to the Kerr one 

, and (c) contains the exact 
A simple way to do that is 



in the test-mass (and test-spin) limit, (b) contains (like -ffroai) the s pin-orbit terms ( |2.' 
spin-spin terms ( 2.54 ) (instead of their "25%" approximation ( 2.51 ) contained in i?reai) 
to define H'^.^^^ = M y^l -|- 2 i^iH'^g — /i)//i with a modified effective Hamiltonian defined as the sum of ( 2.26 ), written 
with the replacement Seff Sq, and of an additional spin-orbit interaction term, A/3'p% with A/3* proportional to 



the difference cr' = SJ: 



On 



Ksix,P,Si,S2) = H^six,p, So) + AHso{x,p,So,cr) , 



where (denoting oq = Sq/M, cos6'o = Sg/ISol) 
AHso{x,p, So,(t) 



At(ao) 



Pi x'' 



(r2-ha2)2-a2 At(ao)sin^ 6*0 



M 



(2.56) 



(2.57) 



with 



— '-'cff ^ '-'0 = 



1/^2 ^.^1^5, 

4 V mi m2 



(2.58) 



The consideration of the new Hamiltonian H'^.^^^ would considerably complicate (even at the qualitative level) the 
discussion of the following section. As we are not sure that this complication really entails a better quantitative 
description of spin effects, when these become important, we shall, in the following, cont ent o urselv es wit h studying 
the consequences of the simpler (though slightly less "accurate") Hamiltonian iJrcai, Eq- ( 2.45 ) with ( 2.46 ). However, 
we mention that it might be useful to consider simultaneously i?rcai and i?rcaiJ ^^'^ to trust their predictions only 
in the cases where they differ only by a slight amount. This gives a useful measure of the domain of validity of the 
present spin-dependent effective-one-body approach. 
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III. DYNAMICS OF TWO SPINNING BLACK HOLES 



A. Equations of motion and exact or approximate first integrals 

In the previous section we have exphcitly constructed an Hamiltonian H^egi{x,p, Si, S2) describing (to some ap- 
proximation) the (conservative part of the) gravitational interaction of two spinning black holes in the center-of-mass 
frame of the binary system. In the present section, we shall describe some consequences of this Hamiltonian. Let 
us sta r t by writing down explicitly the evolution equations for all the dynamical variables. From the basic PB's 
(^l0|)-(|ll^) we get 



= ,3.1) 



(3.3) 



dt ^ 2 ' dSi 
In vectorial notation, the spin evolution equations read (e.g. for the first spin) 

— — = rZi X Si , S2i = „ ^ . (3.5) 
dt obi 

A first consequence of these results is that the magnitudes of the two spins are exactly conserved: 

Si ~ const. , S2 ~ const. (3-6) 

Another general consequence is the exact conservation of the total angular momentum 

J = i + Si + 52 , (3.7) 

where U = e'^^^ pk- Indeed, it is easily checked that J* generates, by Poisson brackets, global rotations of all the 
vectorial dynamical quantities: {J*, V^} = e^^^ for V ~ x, p, Si or 5^2. As the Hamiltonian is a scalar constructed 
out of a;, p. Si and ^2, we have 

^ P ^{J\H,,,,} = 0. (3.8) 
dt 

Therefore 

J* = const, and, in particular, — const. (3-9) 
Evidently, we have also the conservation of the total energy: 

d iJrcal 



dt 



{ffrcal, ^^roal} = =^ iJ,.cal = COnst. (3.10) 



This closes the list of generic first integrals of the evolution. It should be noted that, in general, quantities such 
as or Si ■ S2 are not conserved in time. This means, in particular, that the magnitude of the effective spin, 
— SIq, will not stay constant during the evolution. 

Evidently, in particular situations, more quantities might be approximately conserved. An interesting case is that 
in which the spins are small enough for one to retain only the terms linear in them. In this approximation 

r^2 _ ^(0 = 0)1 

7Ieff(x,p,5i,52) ~iJo(a;,p) + ^ —\ ^ i-5eff, (3.11) 
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where Hq{x,p) is sphericaUy symmetric. 

Let us, more generahy^ assume that H^s, as well as -ffrcai, are spherically symmetric, functions of x and p except 
for a dependence on the combination L ■ ScS- 

-ffroal = H^ca\{r,Pr, , L ■ S cs) , (3-12) 



where Pr = Pi is canonically conjugated to r {{r,pr} = 1). Under the assumption ( 3.12| ) the angular momenta 
evolution equations become (with (LS) = L ■ Scs) 



dSi ^ dHj-cai 
"dT ~ d{LS) 

dS2 _ dH^cal 

"dT ~ d{LS) 

dL _ d Hicai 
~dt ^ c)(LS) 



CTiXxS'i, (3.13) 



CTa-LxS'a, (3.14) 



5'offxZ,. (3.15) 



These evolution equations imply not only (as in the general case) the conservation oi J = L + Si + S2^ and of S\ 
and S^, but also that of: 

= const, L ■ ScS = const. (3.16) 
Note, however, that S'^g is not conserved. Moreover, the radial motion is governed by the equations 

dH,ce,i{r,Pr,L'^,L ■ Scs) ,„-,„^ 
75 , (3.17) 

OPr 

. _ dH,ca.i{r,pr,L'^,L ■ Ses) , . 

Pr — :^ ■ l^o.ioj 

or 

In view of the constancy of = Cq and L ■ ScS = Ci, we see from these equations that the function of r and pr, 
HmdiT^Pr) = -ffroai(^, Pr, C'o: Ci)j dcfincs a rcduccd Hamiltonian describing the radial motion, separately from the 
angular degrees of freedom. In particular, we see (using the fact that Pr enters at least quadratically in -ffrcai) that, 
under our current (approximate) assumption ( ^.12| ), there exists a class of spherical orbits, i.e. of orbits satisfying 

r = const, d H,Mr,Pr ^0, ■ S.s) ^ ^ _ 

or 

Because of the (possibly non-linear) spin-orbit coupling, i.e. the dependence of i?rcai on L ■ Soff, the orbital plane 
of these "spherical" orbits is not fixed in space. But the radial coordinate r being constant, these orbits trace a 
complicated pat h on a sphere (hence the name). These orbits are the analogs, in our two-body problem, and in the 



approximation ( 3.12 ), of similar exact "spherical" orbits for the geodesic motion of test particles in a Kerr spacetime 
p5[ . Their existence (under some approximation) in the two-body problem is interesting for the following reason. 
One expects most black hole binary sources of interest for the LIGO/VIRGO/GEO network to have had the time 
to relax, under radiation reaction, to circular orbits. When the two black holes will get closer, these circular orbits 
will adiabatically shrink until they come close enough for feeling the effect of the spin-orbit coupling (which varies 
oc r~^). In some intermediate domain where the spin-orbit coupling is significant, but couplings quadratic in the 
spins are still small, the initially circular orbit will evolve into an adiabatic sequence of "spherical" orbits of the type 
just discussed. [We are here adding by hand the effect of radiation reaction, treated as an adiabatic perturbation of 
the conservative dynamics discussed in this paper.] These considerations indicate that, in first approximation, the 
total amount of gravitational radiation emitted by coalescing spinning black holes will be determined by the binding 



energy of the Last Stable Spherical Orbit (LSSO), i.e. the last stable solution of Eqs. (3.19), which will satisfy 



For instance, we can assume Eq. (3.11) for //off, but make no further approximation in computing //roai = f{Hcs) 
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d HrcaX I nr2re \ n^^ ^^rcal 

[r,pr = 0, L ,L Scs) = 



di 



di 



(3.20) 



Before studying the energetics of these LSSO's let us mention the existence of other approximate first integrals in 
the dynamics of binary spinning black holes. Let us keep all the terms non-linear in 5'cff , i.e. the full expression of 
Hicai{x,p, Si, S2), but let us try to approximately decouple the orbital motion from the spin degrees of freedom by 
cons idering that the two spin vectors evolve adiabatically (i.e. slowly on the orbital time scale), through Eqs. (t3.3|) , 
{'iA). In this adiabatic-spin approximation, the orbital motion is described by the Hamilton- Jacobi equation (^.171), 



with an adiabatically fixed effective metric . With the definition ( 2.37 ) of the quadratic- in- momenta term Q4{p), 
one can check that, in this approximation, there will exist a two-body analo g of the Carter constant for geodesic 
motion in Kerr ||53|] . Indeed, we have constructed our deformed Kerr metric ( ^.33|) so as to respect its s epara bility 
properties. Let us work in an (adiabatic) Boyer-Lindquist-type coordinate system {t,r,9,ip), as in Eq. (2.33). We 
find that the separability of the effective Hamilton- Jacobi equation yields the following first integrals (of the effective 
Hamiltonian) 



Pt 



(3.21) 



pI 



a Ecff sin 



sm 



+ cos^ = IC=Q+{L,-a E^sf 



(3.22) 



cos 



Q = IC- {L^~aE, 



cS) 



(3.23) 



The last two equations are equivalent to each other, but, d epen di ng on the context, one can be more convenient than 
the other. Let us note the connection of the first integrals ( 0.2l| )-( 3.23|) with the above analy sis o f the first integrals of 
the Hamiltonian depending o nly o n the combination L ■ ScS- The conservation of Lz, Eq. ( 3.2l| ), corresponds to the 
conservation of L ■ SeS, Eq. ( ^.16 ), while the conservation of /C or Q corresponds to the conservation of L^. Indeed, 
if we neglect the terms oc in (3.23) we get 



Q 



'pI 



Li 



1 - sin=^ e 



sm 



(3.24) 



This suggests that, even beyond the adiabatic-spin approximation, the quantities, now defined in an arbitrary frame 
as 



Lz = L-s, Q = L^-{L-sf+a'{n-sf [i? - E^^) , 



(3.25) 



aEcsY) conserved to a good approximation. We are mentionning here 



will be (as well as and K, = Q + [L^ 
these approximate conservation laws because they could be helpful in qualitatively understanding the full two-body 
dynamics. 



B. Spherical orbits and last stable spherical orbits 



We discussed above the existence of spherical orbits under the assumption (or the approximation) that H^^ai depend 
only on the "spin-orbit" combination L ■ 5'cff (as it does at the linear-in-spin level) . More generall y, we have seen that 
if we treat the evolution of the spins as being adiabatic, we have the (approximate) first integrals (3.25). If we use (as 
a heuristic mean of studying the main features of the orbital dynamics) this adiabatic approximation, we can define 
a fam ily of spherical orbits by drawing on the conservation of the quantities (3.25). Indeed, inserting the definitions 
( 3.25 ) into Eq. ( 2.17 ) we get an equation controlling the radial motion: 



At (r) 



A,p2 + 2(4-3i.)i.(GAf)2p4/^2 ^ 
[(r^ + a^) £;,ff - a L.f -{^^r'' + Q+ [L, - a E,s?) ■ 



(3.26) 



The right-hand-side of Eq. ( |3.26 ) defines (when, a? — S'^^/M'^, and Q are considered as adiabatic constants) a 
radial potential whose local minima, in r, determine (adiabatic) spherical orbits. The last stable spherical orbit is 
obtained when this radial potential has an inflection point. More precisely let us define 
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i?(r, E,s, L,, Q) = ((r^ + c?) - a L, f - At{r){^l^ r^ + Q + {L,-a Ef) . 



(3.27) 



The spherical orbits are the solutions of 



R{r, E,s, L„ Q) = Tr R{r, E.s, L,, Q) = . (3.28) 
or 



The solutions of Eq. (3.28) yield a two-parameter family of solutions, along which, for instance, r and E^s are functions 



of Lz and Q. The last stable spherical orbit (LSSO) along such a family of solutions must satisfy the three equations 



R{r,E,s,L,,Q) = — R{r, E,s, L,, Q) = — R{r, E,s, L,, Q) ^ . (3.29) 
or Or'^ 

There is a one-parameter family of LSSO's. For instance, one can take as free parameter the dimensionless ratio Q/L\ 
which is a measure of the maximum angle between the orbital plane and the equatorial plane defined by S'cff. [Note 
that Q = for an orbit in the equatorial plane.] For each value of this angle, and for each value of the effective spin 
parameter a, there will be some LSSO, with particular values of r, i^eff and Lz- 

To study the values of the (effective and real) binding energy, and of the orbital angular momentum along this 
one-parameter family of LSSO's, it is convenient to work with slightly different variables. Let us introduce 

= i, -aSeff, /C= Q + i^ (3.30) 

Let us also work with the radial variable u = \/r and denote 

A{u) ^^^Pi [Al^du) + a' u'] . (3.31) 

We have 

r-'^Rir) = U{u) = {E^s - aLzU^f - A{u){fX^ + ICu^) . (3.32) 
The equation U{u) = (i.e. R{r) = 0) is now solved as 



^^cff = Wa{u, Lz,}C) = aLzu' + ^A{u) {fi^ +ICu^). (3.33) 

The two-parameter family of spherical orbits is now obtained (as functions of the parameters Lz and /C) by solving 
dW/du — 0, while the one-parameter family of LSSO's is obtained by solving dW/du — d^W/du^ — 0. The 
advantage of this formulation is that it exhibits in the simplest way the analogy with the effective radial potential 
discussed in [|lfo| for the pure (3 PN) orbital motion (without spin), namely 



Wa{u,L) = v/A(w)(ai2 + L2 y2)^ (3 34) 

with A{u) = P;^ [A^^{u)]. Apart from the replacement /C, the only two differences between the spinning case 



(( |3.33 )) and the spinless one ( (|3.34 )) is the addition of the spin-orbit energy term + aLzU^, and the additional 



9 2 

a u 



term in the PN expansion of A{u). [Note that A{u) ^ A{u) + a?' v? because the Padeeing is done after the addition 
of a? v? .] We have chosen to parametrize Wa{u) in terms of Lz and K. because it simplifies very much its expression 
and thereby renders more transparent the new physics incorporated in our e ffective one-body approach. The fact 
that Lz depends both on Lz and E^s is not a problem for solving Eq. ( [3.32| ) for E^s- Indeed, we are discussing a 
continuous family of solutions and it is essentially indifferent to parametrize them in terms of L z or L z ■ We could have 



introduced another effective potential W^{r,Lz, Q) by solving R{r, E^s, Lz, Q) = 0, with Eq. ( |3.27D , which would be 
more complicated, but which would describe the same physics. [Note that W^{r) would directly exhibit the correct 
fact that the spin-orbit energy, for given Lz, decreases like r~^, while this fact is hidden in Wa{u) which assumes that 
Lz = Lz~ a EcS is given.] 



C. Binding energy of last stable spherical orbits 



To get a first idea of the physical consequences of our effective one-body description of coalescing spinning black 
holes we have numerically investigated the properties of the one-parameter family of LSSO's. The most important 
quantity we are interested in is the binding energy at the last stable spherical orbit because it is the prime quantity 
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determining the detectability of the GW emitted during the inspiral. We recall that the real, two-body energy is 
related to the effective energy entering the equations of the previous subsection through 



E,^^i = M J 1 + 2iy 



We are mostly interested in the (dimensionless) binding energy per unit total mass, say 




The value of e at the LSSO depends on three dimensionless parameters 



_ . M _ . ^_ a _ \Scs\ 
.4 = 4-^4., a^-^^ 



cose,s^-^. 



(3.35) 



(3.36) 



(3.37) 



Here, the parameter 1/4 (renormalized so that < 1^4 < 1) determines the effect of having comparable masses (1^4 ~ 1) 
rather than a large mass hierarchy (1/4 <C 1). The dependence of e'"^'-' on 1^4 in absence of spins was studied in 
It was found that the ratio e^^'-*/i/4 was essentially linear in 1^4 (even for 1^4 as large as 1, corresponding to the 
equal-mass case) 



oLSO 



-0.014298 1/4(1-^01:^4) • 



(3.38) 



Here, the numerical value 



1.4298% = I 



— 1 ) is one fourth the specific binding LSO energy of a test particle 



in Schwarzschild. The numerical coefficient ci which condenses the effect of resummed PN interactions was found to 
have a value cj^^ ~ 0.048 at 2 PN and c?pN(^^ = 0) ~ 0.168 at 3 PN, and for ojs = 0. [The dependence of cf 
ujs is also roughly linear: cf^'^(a;s) ~ 0.168-1-0.0126^8, at least when — 10<Ci;s<0.] 

We expect that the dependence on V4 of the spin-dependent effects will also be roughly linear (after factorization 
of an overall factor 1/4 which comes from expanding the square root in Eq. ( 3.36 )). In the following we shall generally 
consider (in our numerical investigations) the case ^^4 = 1, and concentrate on the depe ndence on the other parameters. 

Let us clarify the meaning of the parameters a and cos^ls introduced in Eq. (3.37). The quantities S^s, Ez and /C 
entering these definitions are all supposed to be computed at the last stable spherical orbit of an adiabatic sequence of 
spherical orbits (in the sense discussed above). Physically, we have in mind the sequence of inspiralling orbits driven 
by radiation reaction. Technically, we define e(i/4, a, cos^ls) by solving the effective radial potential problem defined 
in the previous subsection. We are aware of the fact that we cannot really attach to cos 0ls the meaning of being 
the cosinus between L and S^s (as the name would suggest), but this is not important. What is important is that 
there is indeed a physical degree of freedom related to the misalignment between Li and 5'off at the LSSO and that 
we measure it by a parameter normalized so that cos^ls — 1 (or —1) when all angular momenta are aligned (in this 
limit the concept of last stable circular equatorial orbit is meaningful and coincides with the cos^ls = 1 (or —1) limit 
of our formal definitions). 

No te that, with the definitions (3.37) and the additional definition 1= i^)^!'^ jGM [i^ the effective radial potential 
(3.33) yields (in dimensionless form, u = GM/r) 



Wa{u, cos Oi^s J) = a cos6'ls^" + J A{u,a^){l + Pu^) 



(3.39) 



This form makes it clear that E^'g^'^ / fi, and therefore 6^"^^°, Eq. ( 3.36| ), will depend primarily on the combination 
("projected value of a") 



a cosf/Ls 



fc ■ ■S'cff 



(3.40) 



where k = L/\L\ (at the LSSO) . For smallish spins, the combination Op is the only one entering the problem linearly. 
As recalled by the notation in ( ^.39 ) the non-projected value of a enters only quadratically in A(u,a^). 
Let us consider more closely the crucial quantity 



1 



(mi 



= 



m2) 
3 



3 7712 

4 mi 



^1 fe • Si 



kSi 



3 mi 

4 m2 



kS2 



^1 ^2 fe • 22 . 



(3.41) 
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In the second form, we have defined Xi = rrii/M, X2 = m2/M {Xi + X2 = 1, X1X2 = i^) and ai = Si/m\, 
0,2 = S2/m\. [We recall that a maximally spinning hole would have |ai| = 1.] 

An important question (for the relevance of the present work) is: what are the plausible values of in the sources 
that will be detected by LIGO/VIRGO/GEO? Present astrophysical ideas about the formation of binary black holes 
suggest neither that the holes be typically maximally spinning, nor that there be any correlation between the 
spin and angular momenta, i.e. between the directions of fe, Si and 82- Not much is known cither about the probable 
value of the mass ratio. To have an idea of the plausible values of Op (which is an algebraic quantity which can 
take positive or negative values) let us consider the random mean square (rms) value of Op under the assumption of 
random, uncorrelated directions fc. Si and S2 (so that (cp) = 0). Let us assume (for simplicity) that mi = m2, i.e. 

— 1, which is the most favourable case because e^^^^^ cx 1^4. We assume also that (Sj^) = (Sj) = (a™ 
given quantity (to be determined by astrophysical models). This yields for a™*' = J (a^) 



a™ = ^ ^ a™^ = 0.357 a™'' . (3.42) 

Even if a™^ = 1 (which would mean that all black holes are maximally spinning) we get a™^ — 0.357. However, 
we find it highly plausible that a™'' will be significantly smaller than 1. For instance, if we optimistically assume a 
uniform distribution of spin kinetic energy between and the maximal value we would get a™^ — and therefore 

a™ — 7/(16 \/3) — 0.253. In view of these arguments, we find plausible that most LIGO/VIRGO binary black hole 
sources will have |ap| <0.3. This consideration is important because we shall see later that for such smallish values 
of flp the simple analytical approach advocated here seems to be quite reliable. However, one should also be able to 
compute physically reliable (or, at least, sufficiently fiexible) templates for fast spinning binary black holes. As we 
said in the Introduction, we think that the EOB approach can be an essential tool for this purpose, in conjunction 
with numerical data, by feeding the (necessarily sparse) numerical data into some multi-parameter version of the EOB 
Hamiltonian. 

This statistical estimate of the plausible value of Op suggests that a typical value of cos^ls — k ■ ^off/ISoffl is 
around ±l/\/3. In our numerical estimates of e'^^^*-' we have used this value, as well as the (implausible) value 
cos^LS = ± 1 corresponding to perfect alignment. As a first step towards exploring the "parametric flexibility" of th e 



EOB approach, we have studied the dependence of e'^^^'-' on the value of the parameter w^, as it appears in Eq. ( ^.2lD . 
We have done numerical simulations for three fiducial values: Ws = = uP'^^ |^l|, u)s — —1987/840 = uo^^ p^,l§||^'], 

and also for uos ^ -\{^ ~ ^^"j = uo*^. Note the numerical values: ujf^ ~ -2.3655, uj^ ~ -9.3439. The original 



2 3 32 ^ 

motivation (when writing this paper, before the completion of the work pT[ | which determined the correct value of w^) 
was to study the sensitivity of our results to the "3PN ambiguity" . We kept it here as an interesting case study of the 
sensitivity of EOB results to modifications of the various coefficients entering the EOB Hamiltonian. The change from 
LOs — Q to LOs = iof^ corresponds to a change of the coefficient 04 = a/i{v)lv from 18.688 to 13.957, i.e. a fractional 
change of —25.32%. The value ojs = has the effect of completely cancelling the 3 FN contribution to the radial 
functions A{u) and A{u). Therefore, choosing uos = gives for the LSSO quantities the same results as the 2 FN 
effective-one-body Hamiltonian Its consideration is useful for exhibiting the difference between the 2 PN-based 
results and the 3 FN-based ones. Our results are displayed in Table ^ and Fig. 1. 

The most important conclusion we wish to draw from these results is that, when a < 0.3 (which, as we argued above, 
covers a large domain of the physically relevant cases) , the binding energy at the LSSO seems to be reliably describable 
by the simple analytical EOB Hamiltonian defined above. Indeed, the differences between: (i) the non-spinning case 
and the spinning ones, and (ii) the 2 FN orbital approximation and the 3 FN one, are all quite moderate (which 
indicates that the effective one-body approach is effective in resumming FN interactions near the LSO). Furthermore, 
the difference between: (iii) the spinning 3 FN case with — — uj^'^^ [|l^, and the same case with — Lof^ 
psf is rather small. This is a testimony of the robustness of the EOB approach. A change of its 3FN coefficient by 
25% does not affect much the physical predictions. This robustness at the 3FN level is indicative of some robustness 
against the addition of higher FN effects. Note also (from table ||) the confirmation that when \ap\ < 0.2, the binding 
energy at the LSSO depends nearly only on the projected effective spin parameter Op = ocos^lSj with a very weak 
dependence on the value of cos 0ls ■ 



^"This value corresponds to taking A = in tj^ = — llA/3 — 1987/840. Here, A denotes the natural ambiguity parameter 
entering the Blanchet-Faye framework. Note that the authors of Refs. p^-[l^ do not claim that A = is a preferred value. 
However, as A is expected to be of order unity we use A = 0, i.e. uis = — 1987/840 as a fiducial deviation from ojs = 0. 
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cos Sls 


= ±l/^/3 






cos Sls 


= ±1 






a 








a 






-0.6 


-1.039 


-0.01319 


6.298 


-0.6 


-0.6 


-0.01150 


7.344 


-0.5 


-0.8660 


-0.01337 


6.220 


-0.5 


-0.5 


-0.01207 


6.989 


-0.4 


-0.6928 


-0.01364 


6.109 


-0.4 


-0.4 


-0.01271 


6.623 


-0.3 


-0.5196 


-0.01405 


5.940 


-0.3 


-0.3 


-0.01345 


6.242 


-0.2 


-0.3464 


-0.01463 


5.700 


-0.2 


-0.2 


-0.01433 


5.841 


-0.1 


-0.1732 


-0.01547 


5.377 


-0.1 


-0.1 


-0.01538 


5.415 


0. 


0. 


-0.01670 


4.954 


0. 


0. 


-0.01670 


4.954 


+0.1 


+0.1732 


-0.01859 


4.391 


+0.1 


+0.1 


-0.01842 


4.439 


+0.2 


+0.3464 


-0.02203 


3.580 


+0.2 


+0.2 


-0.02091 


3.833 


+0.3 


+0.5196 


-0.05146 


1.344 


+0.3 


+0.3 


-0.02529 


3.005 


+0.4 


+0.6928 


-0.1790 


0.7752 


+0.4 
+0.5 
+0.6 


+0.4 
+0.5 
+0.6 


-0.04930 
-0.1048 
-0.1474 


1.538 
1.194 
0.9792 



TABLE II. Binding energies e = (froai/Af) — 1 and (effective Boyer-Lindquist) radii of last stable spherical orbits (LSSO) for 
equal-mass spinning binary systems. The LSSO's depend on two independent parameters: \a\ = \ Scff and cos^LS = L,/y/K. 
(which is, moraUy, the cosine of the angle between the orbital angular momentum and the effective spin). The combined 
parameter Sp = \a\ cos6'ls (projected spin) plays a primary role for moderate spins. The algebraic quantity a is defined as 
+ |a| if cos 9-Ls > and — |a| if cos Sls < 0. All quantities are computed from the 3 PN-level Pade-resummed effective-one-body 
Hamiltonian (p.45[) with los — 0. 
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FIG. 1. Dependence of the binding energy e = {£rca.i/M) — 1 of the LSSO on the effective spin parameter a = Scs/M^ (taken 
with the sign of cos^lsJ- We consider an equal-mass system with a typical misahgnment angle cos6'ls = ± l/\/3. All three 
curves used the Pade-resummed Effective-one-body approach. The lower curves use a 3 PN-level Hamiltonian: the lowest one 
uses LJs ~ = Ljf''^ [10], while the middle one uses coa = —2.365 = ujf^ [17]. The upper curve uses ujs = —9.344 = a;^, which 
is (essentially) equivalent to using a 2 PN-level Hamiltonian [8]. 
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On the other hand, it must be admitted that when, say, a> 0.4 the differences between the three cases (i), (ii), 
(iii) become so large, and the radius of the LSSO becomes so small, that the present spin-dependent EOB predictions 
cannot be quantitatively trusted. [However, as discussed in more detail below, we think that they remain qualitatively 
correct.] If the orbital dynamics were well described by the 2 PN-level orbital EOB metric (i.e. if ajg had been near 
—9; see upper curve in Fig. |l|), the binding energy, even in such extreme cases, would differ only moderately from 
the non-spinning case, and we could trust the EOB-plus-spin predictions. However, as uj^ is zero the 3 PN EOB + 
spin predictions become, for cos6'ls — 1 and a>0.4, very different from the 2 PN ones and quite sensitive to the 
numerical values of the expansion coefficients entering the EOB potentials. Let us note, however, that in all cases 
(even the most extremely spinning ones) if the spin parameter Op is negative (i.e. if the effective spin vector, whatever 
be its magnitude, has a negative projection on the orbital angular momentum) the EOB predictions become extremely 
reliable because all the differences between the cases (i), (ii), (iii) become quite small. 

All these results are easy to interpret physically. This can be seen from the basic equations of the EOB approach 
which simplify so much the description of the physical interactions by representing them as slightly deformed versions 
of the well-known gravitational physics of test particles in Schwarz schild or Kerr geometries. Indeed, the basic equation 
of the EOB approach determining the binding of the LSO is Eq. ( 3.39| ) which differs from its wcU-knownp] Kerr limit 
(i.e. v 0) only by the change 

AK{u,a'^) = l- 2u + a^u'^ ^ A(u,a^) = [1 -2u + c?u^ + 2vu^ + ai{v)u^] . (3.43) 

The crucial point (which is, finally, the most important new information obtained by the 2 PN and 3 PN orbital 
calculations) is that the 2 PN and 3 PN additional terms to the radial function A™ q(u) have both positive coefficients. 
This means that, even before the addition of the effect of spin (which leads to a -|-a^u^ additional term in Ax{u), 
corresponding to the famous + term in Ax{r) — — 2Mr + a^) the main effect of non-linear orbital interactions 
for comparable masses is "repulsive", i.e. correspond to a partial screening of the basic Schwarzschild attractive 
term 1 — 2m = 1 — 2GM/c^r by the addition of repulsive terms oc +i'/r^ and +i'/r^. Now, paradoxically, the 
addition of a repulsive term leads to a more tightly bound LSSO because the less attractive, but still attractively, 
radial function A{u,a'^ = 0) will be able to "hold" a particle in spherical orbit down to a lower orbit. In other words, 
when a radial potential becomes less attractive, its LSSO gets closer to the horizon, and the binding energy of the 
LSSO becomes more negative. This being said, one understands immediately the additional effects due to the spin 



interaction. There are basically two such effects: (a) a linear "spin-orbit" effect linked to the +apiu'^ term in (3.3£) 
(with Op = acos^Ls): and (b) a non-linear spin-quadratic modification of the metric coefficient, i.e. the additional 
+ 0^ term in A{u,a^) {or in Ak{u) = l — 2u + a^iP). The crucial points are that: (1) when < 0, i.e. cos^ls < 
(coarse antialignment of angular momenta) the dominant linear spin-orbit coupling is attractive and therefore pushes 
the LSSO upwards, towards a less bound orbit, while, (2) when Op > 0, i.e. cos6'ls > (coarse alignment of angular 
momenta) both the linear spin-orbit coupling +ap£u'^ and the spin-quadratic additional term -I- v? are repulsive 
and tend to draw the LSSO downwards, i.e. closer to the horizon, in a more bound orbit. Therefore we see that, when 
Op > 0, all the new effects (the z/-dependent non- linear orbital interactions and the spin effects) tend in the same 
direction: towards a closer, more bound orbit. As the existence of a LSSO is due to a delicate balance between the 
attractive gravitational effects and the usual repulsive ("centrifugal") effect of the orbital angular momentum (i.e. the 
term + Pu^(x+ /r"^ in (3.39)), when several attractive effects combine their action, they start having a large effect 



on the binding of the LSSO. This is well-known to be the case for circular, equatorial, corotating (Sp — ~\-'a) orbits of a 
test particle in Kerr, which feature, in the case of an extreme Kerr (a = 1) an LSO at r = 1, with /i-fractional binding 
(i?off — Ai)/M = l/VS — 1 = —0.42265 (corresponding to e ~ v {Ecu — /^)//^ — —0.10566 1^4). It is also well-known 
that, again for extreme Kerr, a counterrotating (cp = — a) circular, equatorial orbit in extreme Kerr has an LSO at 
r = 9, with /i-fractional binding (i?eff — m)/m = 5/(3-^/3) ^ 1 = —0.037750 (corresponding to e ~ —0.0094374 z/4). 
What is less well-known is that the extreme binding of the circular, equatorial, corotating LSO around an extreme 
Kerr is not representative of the binding of typical LSSO's around typical (or even extreme) Kerr holes. Indeed, when 
cos^LS 7^ i 1 (i-e. in more invariant language, when Q 7^ 0) and when a < 1, the LSSO is, in general, moderately 
perturbed away from the Schwarzschild value tlso 6 and its binding is correspondingly moderately different from 



Actually, as far as we know, the Kerr limit of Eq. (15.39) has never been written down before. Usual treatments use the 
more complicated effective radial potential Wa{T, Lz,Q). Anyway, the physics is the same, but it is more cleanly presented in 

(ip- 

^ Remember that we Fade resum A{u) and A{u) to ensure that these functions qualitatively behave like 1 — 2 u or 1 — 2 
(for < 1), i.e. (generically) have a simple zero near r = u~^ = 1-1- V 1 — 2^ , which means that the effective metric becomes 
"infinitely attractive" at some deformed horizon. 
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its Schwarzschild limit eg^j^ ~ —0.014298 ;/4. The present work has shown that the location of the LSSO's for binary 



spinning holes can be rather simply obtained, in the EOB approach, by balancing in the specific way of Eq. (3.39) the 
centrifugal effect of the orbital angular momentum against the overall attractive effect of gravity, but with the critical 
addition of the 2 PN and 3 PN repulsive terms, of the spin-quadratic repulsive term, and of the indefinite-spin effect 
of the spin-orbit interaction. 

D. Expected spin of the hole formed by the coalescence of two spinning holes 

The last topic we wish to discuss concerns the expected result of the coalescence of two holes. In particular, we 
are interested in estimating the maxim al sp in that the final hole, resulting from the coalescence of two spinning 
holes, might have. It was estimated in 11| (by using the EOB approach) that the coalescence of two non-spinning 



holes of the same mass mi = m2 = M/2 leads (after taking into account the effect of gravitational radiation on the 
orbital evolution and on the loss of energy and angular momentum) to the formation of a rotating black hole of mass 
Mbh — (1 — erd) 0.976 M and spin parameter obh — 0.80 [we have included a factor 1 — e^d in Mqh to take into 
account the energy loss during the ring-down. Ref. jl^ found e^d — 0.7%]. The fractional energy 0.976 — 1 = —0.024 
roughly corresponds to the (adiabatically estimated) LSO binding energy (—0.015 in the 2 PN-based estimate of |^) 
minus the energy per unit mass radiated during the plunge (~ —0.007 ||ll|). We shall leave to future work a similar 
estimate, for the 3 PN-plus-spin case, of the amount of energy emitted in GW. We wish here to focus on the issue 
of the spin of the final hole. The above value obh — 0.80 is rather close to the maximal value flmax — 1 ^^nd there 
arises the question of whether an EOB treatment of the coalescence of two spinning holes might not formally predict 
a final value of obh larger than one! By "EOB treatment" we mean here a completed version of the EOB approach 
(as in at the 2 PN, non-spinning level) obtained by: (i) adding a resummed radiation force to the "conservative" 
EOB dynamics, and (ii) pushing the calculation of the EOB evolution down to its point of unreliability (near the 
last unstable orbit) where it is matched to a perturbed-single-black-hole description. A zeroth approximation to this 
completed EOB approach is the one we study in this paper: an adiabatic sequence of solutions of the conservative 
dynamics, terminated at the LSO. In this approach one entirely neglects the losses of energy and angular momentum 
during the plunge phase following the crossing of the LSO. The numbers recalled above show that the energy loss 
during the plunge (and the ring-down) is not negligible compared to the binding energy at the LSO. However, for the 
present question this is not a problem. What is important is that the angular momentum loss during plunge is a very 
small fraction (a percent or so) of the angular momentum at the LSO, and that the final mass of the black hole is 
nearly equal to M = mi + m2 . This leads us to the following zeroth order estimate of the spin parameter of the final 
hole: 

. |J|LSSO |j|LSSO 

«BH-^^^-^^. (3.44) 

In view of the exact conservation of J in our conservative EOB (real) dynamics, it is clear that it is | Jj''^^^ which is 
a good measure of the total angular momentum of the final spacetime, i.e. of the final black hole. 



We are facing here a potential consistency problem of this simple-minded EOB treatment: when computing ( 3.44 ) 
for spinning configurations does one always get obh < 1? One might worry that, starting with a value of obh — 0.80 
for non-spinning holes, the addition of large spins on the holes might quickly exceed the extremal limit. It is plausible 
that the most dangerous situation is the "aligned c ase" , where all the angular momenta, L, Si and 5'2 are parallel 



(or antiparallel) . In this case the numerator of Eq. ( 3.44 ) reads 

jLSSO ^ ^LSSO + S1+S2, (3.45) 

while the spin parameter of the effective metric reads 

« = «p - - (^1 + r) + (^^' + r) ■ ^^-^^^ 

Here, we consider S'l, S2 and ai = Si/m\, 0,2 = 52/m2 as algebraic numbers (positive or negative). This allows us to 
investigate also the case where the spins might be antiparallel to k. For simplicity, we shall only study the symmetric 
case where mi = m2 and S'l = 6*2. For this case 

7_ 

a = flp = - ai , (3.47) 

and 
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a 




FIG. 2. Approxima.te prediction for the spin parameter Sbh — \ J\^^^'^ /M^ of the blaclc hole formed by the coalescence 
of two identical spinning holes (with spins parallel or antiparallel to the orbital angular momentum). The horizontal axis is 
the effective spin parameter a — |ai = |o2. The three curves correspond to the three cases plotted in Fig. Note the 
prediction (robust under changing the 3 PN contribution to the effective potential by 25 %) that the final spin parameter is 
always sub-extremal, and reaches a maximum Sbh — 0.87 for a ~ + 0.3. 



tLSSO 1 1 1 4 

~ ~ 4 ^ 2 4^ + 7 (.■^■4^^) 

where the dimensionless orbital angular momentum = L^/^M is related to the dimensionless quantity (when 
cos 0LS = 1)1= ^/K/iJiM = L,/nM appearing m ( |3.39D through 



L,=i + a — . (3.49) 
It is interesting to note that, even in the case where both holes are extreme (oi =02 = 1) the maximu m va lue of 



the effective spin parameter is amax = | < 1- We have numerically investigated the quantity qbh, Eq- ( 3.48 ), as a 
function of the effective a = ap. The result is plotted in Fig. ^ for different values of the 3 PN parameter uJs- 

We see that the final spin parameter reaches a maximum for a positive value of Op, i.e. for parallel (rather than 
antiparallel) spins. For the correct 3PN value tOs = the maximum value of Sbh is comfortably below 1: namely, 
jj-jmix ^ 0.87, reached for ~ + 0.3. This is not much larger than the value qbh — 0.82 obtained for Up = 0. We find 
that this is a nice sign of the consistency of the EOB approach. This consistency was not a priori evident. In fact for 
UJs < —9 one gets a maximum value of obh slightly larger than 1. In particular, note that the 2 PN treatment of the 
orbital dynamics (obtained for = ~ —9.3439; upper curve in Fig. |^) formally leads to problematic over-extreme 
values of obh- This may be interpreted as a confirmation of the need of "repulsive" 3 PN effects (i.e. + 9 ^ 1). It 
is (a posteriori!) easy to understand physically why (Ibhi after reaching a maximum, then decreases when one adds 
more spin on the two black holes. Indeed, there is here a competition between two effects: adding spin on the holes 
(i.e. increasing Op), on the one hand directly contributes to augmenting obh through the second term of the RHS of 



Eq. ( |3.48 ), but, on the other hand, indirectly contributes to reducing the total J^^sso reducing L^^^*-* (indeed, as 



we explained above, positive spin leads to an LSO orbit closer to the horizon, and therefore with less orbital angular 
momentum). The first effect wins for smallish spins, while the second (more non-linear) effect wins for larger spins. 



IV. CONCLUSIONS 



We started by recalling the need of techniques for accelerating the convergence of the post-Newtonian (PN) ex- 
pansions in the last stages of the inspiral of binary systems. We summarized the evidence (Table |) showing the 
remarkable convergence properties of the best current resummation technique: the effective one-body (EOB) ap- 
proach of Refs. I^,]!^. We showed how to generalize the EOB approach to the case of two spinning black holes 
with comparable masses = fi/M ~ 1/4). As a first step towards computing the spin-dependent EOB Hamil- 
tonian we constructed an effective metric, which can be viewed either as a i^-deformation of the Kerr metric or 
as a spin-deformation of the z^-deformed effective metric. The effective spin entering this deformed Kerr metric is 
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M a = SoS = + I ^) 'S'l + ^1 + I 82- The introduction of this effective a aUows one to combine in a simple 

manner all (PN leading) spin-orbit coupling effects, and most of the spin-spin ones, with the rather complex but 
important 3 PN effects, which have been incorporated only recently in the EOB app roach llG]. 



We have also constructed a more complicated modified effective Hamiltonian, Eq. ( 2.56 ), which separately depends 

on two (effective) spin vectors, Mao = So = (l + ^) S, + (l + ^) S2, and ^ = S,«-So = -i 5i + ^ ^2) , 

and which allows a (hopefully) more accurate representation of spin-spin effects. We recommend the simultaneous 
consideration of Hcs and H^g to determine the domain of trustability of the presently constructed spin-dependent 
EOB Hamiltonian. Namely, when ifroai — f{Hcs) and i?^cai — fi-^cs) ^^ad to numerically very similar evolutions, 
one is entitled to trusting them both; while a significant difference in their predictions signals a breakdown of the 
trustability of the simple EOB Hamiltonian proposed here. 

The present paper has only investigated a few aspects of the physics predicted by our spin-generalized EOB ap- 
proach. In particular, as a first cut toward understanding the relevance of our construction for gravitational wave 
(GW) observations we have discussed the approximate existence of "spherical orbits" (orbits with fixed radial coordi- 
nate, as in the Kerr metric) and we studied the binding energy of the last stable spherical orbits (LSSO). A message 
of this study is that, for most physically relevant cases (in the parameter space where one randomly varies all angles 
and all spin values), the results are only weakly dependent on the exact numerical values of the 3PN coefficients. 
Moreover, they exhibit moderate deviations from the non-spinning case (see Fig. 1 and Table II) . To give a numerical 
fiavor of the effects of spin we note that, when the projected spin parameter ap = k ■ a/M, Eq. ( ^.41 ), is smaller than 



about -t-0.2, its effect on the fractional binding energy (e = (ii'roai — M)/M) of the LSSO is, approximately, 

100 e^^^° ~ -1.43 i^4(l + 0.168 1^4,) ~ 0.806 1^4(1 + 0.888 ap , (4.1) 

where 1/4 = 4;^ = 4mi77i2/(mi -I- 7712)^ < 1. As in most cases (random angles, random spin-kinetic energies) it is 
plausible that |ap|<a™ ~ 0.25, we expect that spin effects will only modify the energy emitted as gravitational 
waves up to the LSSO by less than about 0.6% 1/4 M. 

Such an increase, though modest, is still a significant fractional modification of the corresponding energy loss 
predicted for non-spinning systems (cq = — 1.67%M for ^^4 = 1). In fact, this effect might cause an important bias in 
the first observations. If the intrinsic spins of the holes can (at all) take large values, the highest signal-to-noise-ratio 
events in the first years of LIGO observations might select binary systems with rather large and rather aligned spins. 
It is therefore important to include spin effects in the data analysis of coalescing black holes. We have argued that, 
in most cases, the simple-minded generalized EOB approach presented here should be a reliable analytical tool for 
describing the dynamics of two spinning holes and for computing a catalogue of gravitational waveforms, to be used 
as matched filters in the detection of GW's. However, it must be admitted that, in the cases where the effective 
spin vector is coarsely (positively) aligned with the orbital angular momentum, and where the spins are so large that 
a > 0.4, the predictions from the above-introduced EOB Hamiltonian start predicting LSSO radii so near the "effective 
horizon" where At(r) = 0, that they cannot be quantitatively relied upon. [T hough I would still argue that they can 
be qualitatively trusted, in view of the simple physics they use; see subsection III C| .] We give some examples of that 



m Table 0. In such cases the EOB approach does predict much larger energy losses, possibly larger than 10% M. 
In these cases, the uncertainty in the waveform may be so large that one may need the type of non-linear filtering 
search algorithm advocated in Ref. We wish, however, to emphasize the differences between our treatment and 
conclusions and those of Flanagan and Hughes. These authors defined the "merger" phase as (essentially) what comes 
after the binary system crosses the non-spinning LSO (around 6 GM), and they assumed that the signal from the 
"merger" phase can only be obtained from numerical relativity. Moreover, they optimistically assumed that (in all 
cases) 10% M are emitted in GW energy during the merger phase, and 3% during the subsequent ring-down phase. 
By contrast, our treatment is based on the idea that a suitable resummed version of the PN-expanded dynamics, 
namely the EOB-plus-Pade approach, can, in most cases, give an analytical handle on the computation of the inspiral 
signal down to the spin-modified LSSO (and even during the subsequent plung e, as discussed for the spinless case in 



). We have argued in several ways that the simple EOB Hamiltonian (2.45) gives reliable answers in most cases, 
and allows one to analytically control the possible amplification (or deamplification, when (L < 0) in GW energy loss 
due to spin effects. Moreover, it is only in rather extreme cases that we could agree with Q in predicting > 10% M 
energy losses. In most other cases, we think that the EOB method provides a reliable basis for computing families 
of wave forms that will be useful templates for the detection of GW's. Another difference with |^ is that we have 
argued, on the basis of definite computations, that the spin of the final hole will never become nearly extremal (even 
if the initial spins are extremal). This is important for the data analysis of the ring-down signal, because the decay 
time of the least damped quasi- normal-mode starts becoming large only for near extremal holes. 

Let us emphasize that the present work is only a first step toward an improved analytical understanding of the 
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last stages of inspiral motion of two spinning compact object|^. The explicit spin-dependent Hamiltonian (2.45) (or, 
better, H^.^^^ defined at the end of Section II) has only taken into account the leading effects (in a PN expansion) of 
the spin-dependent interactions. More work is needed to analytically determine more accurate versions of the EOB 
Hamiltonian. In particular, it would be interesting to explicitly derive the next-to-leading ( 1 PN) corrections to the 
EOB spin-orbit (and spin-spin) interactions. Furthermore, we have only provided a resunimation of the conservative 
part of the dynamics. There remains the important complementary task of resumming the radiation reaction part. 
This was done in |^ , using previous results of 0] , only for the spinless case. 

Once this is done, we expect, as in our previous study that the presence of a LSSO along the sequence of 
adiabatic orbits will be blurred and will be replaced by a continuous transition between inspiral and plunge. There 
remains also the task of studying the effects of spin-dependent interactions on the gravitational waveform emitted 
during the last stages of inspiral and during the plunge (that we have not considered here). In other words, one 
needs to redo, by combining the EOB approach with resummed versions of radiation reaction, the studies, valid 
far from the LSO, which were based on straightforward PN-expanded results [Q, Note that our result above 

about the primary importance of the single parameter Op, combined with the understanding that the number of 
"useful" cycles in the GW signal for massive binaries is rather small, suggests that a rather small number of "spinning 
ternplates" will be really needed in a matched filter data analysis. On the other hand, we recall that it was found 
in H that the plunge signal (but not the ring-down one, for stellar mass holes) plays a significant role in the data 
analysis. 

It will also be interesting to see, within the EOB approach, the extent to which the non-linear spin-dependent 
interactions might, as has been recently suggested |Q, lead to a chaotic dynamical evolution. We a priori suspect 
that two factors will diminish the significance of such chaotic evolutions: (1) they occur only in an improbably small 
region of phase space (involving, in particular, large spins), and (2) their effect on the crucial GW phasing is rather 
small. 

It would be very useful to have independent means of testing the accuracy of the EOB approach. At this stage we 
see only three ways of doing that (beyond the performance of more internal checks of the robustness of the approach): 
(i) an analytical calculation of the 4 PN interaction Hamiltonian, (ii) a comparison between numerical computations 
and the EOB results, and/or (iii) a comparison between the EOB predictions and the forthcoming GW observations, 
(i) would be important for assessing the convergence of the PN-resummed EOB Hamiltonian. In view of the extreme 
difficulties involved in the 3 PN calculations ||l^-|l8|] it would seem hopeless to even mention the 4 PN level. But 
in fact, the EOB approach itself suggests that the current methods used in PN calculations are highly inefficient, 
and unnecessarily complicated. Indeed, as emphasized in the final, gauge-invariant content of the 3 PN result 
is contained in only three quantities 04, 63 and Z3, and only one of them, 04(1/), is really important for determining 
the dynamics of inspiralling quasi-circular binaries. If one could invent a new approximation scheme which computes 
directly 04 (at 3 PN), it might be possible to compute its 4 PN counterpart, 05(1^). (ii) is not yet possible (at least as 
a test of the EOB Hamiltonian) because numerical computations use as initial data geometrical configurations that do 
not take into account most of the crucial physics incorporated in PN calculations. Current numerical computations use 
somewhat ad hoc "binary-black-hole-like" data, often of the restricted spatially conformally flat type, without trying 
to match their initial data to the near LSO configurations predicted by (resummed) analytical approaches. On the 
other hand, let us stress that the value of the radial PN-expanded potential A{u) = 1 — 2u + a3{i') + 04(1') u"^ + ■ ■ ■ 
crucially depends on the non-linear gravitational interactions linked to the hfj^ part of the spatial metric, i.e. to its non 

conformally flat part, and also to the non-linear interactions linked to the tt^j, part of the gravitational field momenta. 
For instance already at the 2 PN level, the truncation of the Einsteinian prediction for the two-body problem (driven 
into a close orbit by a long past interaction involving retarded GW interactions) corresponding to artificially assuming 
a conformally flat spatial metric changes the physically correct value a§^^(i^) = 2v into a§°"^ ''''*(z/) = j (18 — 5^) v 
pO|. [It also slightly changes the energy map /.] For equal- mass systems, this corresponds to multiplying the positive 
a^v) by a factor 4-2.09375. As we discussed above, this (artificial) increase of the "repulsive" character of the non- 
linear gravitational interactions tends to artificially increase the binding of the LSO. As the 2PN coefficient is anyway 
too small to have a large impact on the LSO characteristics, this 2PN change does not, by itself, change much the 
LSO energy . However, if this tendency to increase the "repulsive" character of the PN expansion (caused by the 
neglect of the hfj^-, and Tr^^-dependent interactions) persists at the (numerically more important) 3PN level, this 
might explain the current discrepancy between analytical and numerical estimates of LSO characteristics. In fact. 



Though, in most of the paper we only spoke of binary black holes, it should be clarified that our EOB Hamiltonian also 
applies to binary spinning neutron stars or to spinning neutron-star-black-hole systems, at least down to the stage where the 
quadrupole deformation of the neutron star becomes significant. 
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we note that the initial data taken by a recent attempt |59| at fulfiUing the proposal of ||^ to start a full numerical 
calculation only at the moment where it is really needed, i.e. after crossing the LSO, uses LSO initial data |6^] 
with a binding energy clso = E/{mi + 7/12) — 1 ~ —2.3% which is 38% larger than the value e£gQ ~ —1.67% 
obtained at 3 PN (with tOg — 0) by analytical estimates. Similarly, the LSO orbital period of the initial data of |5^] 
is Tlso - 35 (mi + 7722) (6O), which is twice smaller than the 3 PN estimate T^g™ - 71.2 (toi + TO2) 0! These 
discrepancies between state-of-the-art numerical LSO initial data and state-of-the-art analytical estimates of LSO data 
are significantly larger than the natural "theoretical error bar" on the (resummed) analytical estimates (derived, say, 
by comparing 2PN estimates to 3PN ones). [See, however, the new numerical approach of ||6^ whose LSO data agree 
well with the 3 PN EOB estimates ||6^ .] In our opinion, this makes it urgent for the numerical relativity community 
to develop ways of constructing initial data that correctly incorporates the crucial non-linear physics (linked to the 
hjj'" and Tr^y parts of the metric) which is taken into account in PN calculations. If a significant discrepancy remains 
after this is done, one will be entitled to blame the lack of convergence of the EOB-resummed PN calculations. If 
one finds agreement, this will be a confirmation of the claim made here that the Pade-improved EOB is a reliable 
description of the last orbits before coalescence. Once one succeeds in matching analytical and numerical results for 
non-spinning black holes, it will be very interesting to use numerical data on fast-spinning black holes to refine the 
EOB Hamiltonian by fitting the values of the extra parameters which can be introduced in the EOB Hamiltonian to 
represent higher PN effects. In the long term we expect that such a complementarity between numerical and analytical 
tools will be needed for defining a sufficiently dense set of GW templates. [In view of the large dimensionality of the 
parameter space of the two spinning hole system, it seems hopeless to use only numerical techniques to define a dense 
network of templates.] 

Finally, even if no decisive progress is made on (i) or (ii) before the first sources are detected, there remains the 
possibility that the first observations might confirm the soundness of (or suggest specific modifications of) the EOB- 
based waveforms, and thereby facilitate further detections by narrowing the bank of templates. For instance, one 
might include a 4PN contribution +a5{i')u^ to A(u), as a free parameter in constructing a bank of templates, and 
wait until LIGO /VIRGO / GEO get high signal-to-noise-ratio observations of massive coalescing binaries to determine 
its numerical value. 
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